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Abstract 

In this paper we study a continuum version of the Potts model, where particles are 
Q>s^ points in R'^, d > 2, with a spin which may take S > 3 possible values. Particles with 

^-H different spins repel each other via a Kac pair potential of range 7"^, 7 > 0. In mean 

field, for any inverse temperature (3 there is a value of the chemical potential at which 
S + 1 distinct phases coexist. We introduce a restricted ensemble for each mean field 
pure phase which is defined so that the empirical particles densities are close to the mean 
field values. Then, in the spirit of the Dobrushin Shlosman theory, [9], we prove that 
while the Dobrushin high-temperatures uniqueness condition does not hold, yet a finite 
size condition is verified for 7 small enough which implies uniqueness and exponential 
^ decay of correlations. In a second paper, [8], we will use such a result to implement the 

1—1 Pirogov-Sinai scheme proving coexistence of 5 + 1 extremal DLR measures. 
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1 Introduction 

In this paper we consider a continuum version of the classical Potts model, namely a system 
of point particles in where each particle has a spin sG{1,..,5},S'>1, and particles with 
different spins repel each other, this being the only interaction present. When S = 2 this is a 
simple version of the famous Widoni-Rowlinson model which has been the first system where 
phase transitions in the continuum have been rigorously proved, [15], and for S > 2 and at 
very low temperature, a phase coexistence between the S symmetric phases for continuum 
Potts models was established in [11]. 

The mean field version of the continuum Potts model has been recently studied in [10]. The 
phase diagram has an interesting structure. In the A)-plane, /3 the inverse temperature, 
A the chemical potential, there is a critical curve, see Figure 1, above which (i.e. A "large"), 
there is segregation, namely there are S pure phases, each one characterized by having "a 
most populated species" (of particles with same spin). Instead, below the critical curve there 
is only one phase, the disordered one where the spin densities are all equal. The behavior on 
the critical curve depends on S. li S = 2 there is only the disordered phase while if 5 > 2 
there is coexistence, namely there are S" + 1 phases, the "ordered phases" where there is a 
spin density larger than all the others and the disordered phase as well. 
An analogous phenomenon occurs in the mean field lattice Potts model where at a critical 
temperature there is a first order phase transition with coexistence of 5+1 phases if S > 2, but 
in the continuum there is an extra phenomenon occurring at the transition, namely the total 
particles density has a strictly positive jump when going from the disordered to an ordered 
phase. This can be seen as an example of interplay between magnetic and elastic properties 
and interpreted as a magneto-striction effect, as the appearance of a net magnetization is 
accompanied by an increase of density and thus a decrease of inter-particles distances. 
Our purpose is to prove that the above picture remains valid if mean field is replaced by a 
finite range interaction. Let q = (..., rj, s^, ...), i = l,..,n, Vi G M'^, Si G {1,..,5}, a finite 
configuration of particles. We suppose that their energy is 




(1.1) 



where 




(1.2) 
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Jj{r,r') = 7*^7(0, 7(r' — r)), 7 > a Kac scaling parameter, J(0,r) a smooth probability 
kernel supported by \r\ < 1/2. (Observe that Hx{q) is independent of the particles labeling). 
To motivate the above choice recall that the mean field energy density (mean field energy 
over volume) is 

s^s' s 

where p{s) is the density of particles with spin s. Then 

Hx{q)= i ex{pq,r{-))^ PqAs) = Y'^si=sJ-i{r,ri) (1.4) 

i 

Thus H\{q) is the integral of the mean field free energy density, where the latter is computed 
using the empirical averages pq^r{s)- If 7 is small one may think that (1.1) "simulates mean 
field" . Indeed we will prove in [8] that 

Theorem 1.1. For any d > 2, S > 2 and /3 > there is 7* > such that for any 7 < 7* 
there exist Xp^^ and S + 1 mutually distinct, extremal DLR measures at (/?, A^^-y). 



To keep the statement simple we have not reported all the information we have on the 
structure of the DLR measures referring to [8] for the full result. In particular we know 
that the particles densities are close to their mean field values (for 7 small). The proof 
of Theorem 1.1 follows the Pirogov-Sinai strategy which is based on the introduction of 
"restricted ensembles" where the original phase space of the system is restricted by constraints 
which impose local closeness to one of the putative pure phases, in our case local closeness 
of empirical averages to the mean field values in a pure phase. We need a full control of 
such "restricted ensembles" and then a general machinery applies giving the desired phase 
transition. As a difference with the classical Pirogov-Sinai theory, here the small parameter 
is the inverse interaction range 7 instead of the temperature, as we are "perturbing" mean 
field instead of the ground states, see for instance the LMP model, [13], where these ideas 
have been apphed to prove phase transitions for particles systems in the continuum with Kac 
potentials. 

In the typical applications of Pirogov-Sinai, restricted ensembles are studied using cluster 
expansion which yields a complete analyticity (in the Dobrushin-Shlosman sense, [9]) charac- 
terization of the system. Namely constraining the system into a restricted ensemble raises the 
effective temperature and the state enjoys the characteristics of high temperature systems. 
An analogous effect has been found in the Ising model with Kac potentials, [7], [5], and in the 
LMP model, in both the high-temperatures Dobrushin uniqueness condition has been proved 
to hold. This is a "finite size" condition, and the Dobrushin uniqueness theorem states that 
if such a condition is verified, then there is a unique DLR state. The importance of the result 
is that the condition involves only the analysis of the system in a finite box: loosely speaking 
it is a contraction property which states that compared with the variations of the boundary 
conditions, the Gibbs measure has strictly smaller changes, all this being quantified using the 
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Wasserstein distance. Dobrushin's high temperatures means that the size of the box [where 
the conditional measures are compared] can be chosen small (a single spin in the Ising case) 
or a small cube in LMP so that there is no self interaction in Ising or a negligible interaction 
among particles of the box (in LMP) and the main part of the energy is due to the interaction 
with the boundary conditions. The measure and its variations are then quite explicit and it 
is possible to check the validity of the above contraction property. 

As explained by Dobrushin and Shlosman, one expects that when lowering the temperature 
the above high temperature property eventually fails, the point however being that it could be 
regained if we look at systems still in a finite box but with a larger size, eventually divergent 
as approaching the critical temperature. The problem is that if the finite size condition 
involves a large box then self interactions are important and it is difficult to check whether 
the condition is verified. 

While it is generally believed that the above picture is correct, there are however not many 
examples where it has been rigorously established. Unlike Ising with Kac potentials and 
LMP, in an interval of values of the temperature, where the high temperature Dobrushin 
condition is valid in restricted ensembles, in the continuum Potts model we are considering 
there is numerical evidence (at least) that it is not verified. We will prove here that a finite 
size condition (involving some large boxes where self interaction is important) is verified 
in our restricted ensembles and then prove using the disagreement percolation techniques 
introduced in [2] , [3] , that our finite size condition implies uniqueness and exponential decay 
of correlations and all the properties needed to implement Pirogov-Sinai, a task accomplished 
in [8]. 



Part I 

Model and main results 



2 Mean field 

The "multi-canonical" mean field free energy is 

F"''{p) = ^J2p^P^'-^IY.P^^^''^P^-'^^^ P={pi,-,Ps}eRl (2.1) 

s^s' ^ s 

where ps represents the density of particles with spin s and /3 the inverse temperature, to 
underline dependence on P we may add it as a subscript. The "canonical" mean field free 
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energy is instead 

r'{x) = inf (F-f (p); J2ps = ^}^ ^ > (2-2) 

s 

and the mean field free energy CEf^^{x) is the convex envelope of F^^^{p), 

and CEf^^{x), A G M the chemical potential, are defined by adding the term —Ax, where in 

the case of F^^{p), x = ^ ps- 

s 

Observe that for any a > 0, 

F^i{p) = a-'F^/^^,,{ap), A = a-^' - ^ (2.3) 

so that if the graph of CEfp^j^{x) has a horizontal segment, then for any /?', CEfp}^y{x) has 
also a horizontal segment when A' = oA + /3~^aloga, a = /?//?', which reduces the analysis of 
phase transitions to a single temperature, object of the following considerations. 
As shown in [12] (see the proof of Theorem A.l therein), the variational problem (2.2) is 
actually reduced to a two-dimensional problem because: 

Lemma 2.1. 

r'{x) = inf (F-f (p); ^ p{s) = x; p, > p2 = ■ ■ ■ = ps} (2.4) 



The analysis of (2.4) yields: 



Theorem 2.2. Let S > 2 and P > 0. Then there are < X- < x+ such that CEf^^{x) 
coincides with /™^(a;) in the complement o/(x_, a;+) and it is a straight line in As a 

consequence there is \p such that CEf^j^^{x) has the whole interval [x-,x+] as minimizers, 
it is strictly convex in the complement and D"^ fj^j^^{x±) > 0. 



By using the scaling property (2.3) we then obtain the phase diagram in Figure 1. 
We will next discuss the structure of the minimizers of Fp^j^^{p). 

Theorem 2.3. Let S > 2, (3 > and Xp as in Theorem 2.2. Then F^l^{p) has S + 1 
minimizers denoted by p^^\ fc = 1, .., 5 + 1. For k < S, p^^^ > pi''\s / k and p^J'^ = p^^'' for 
all s, s' not equal to k. Lnstead pi^^^^ = Pi^'^^^ for all s and 

s s 
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Figure 1: Phase Diagram of the Mean field Potts gas 

Finally for any k the Hessian matrix L^'^^ := D'^F"^j^^(p^^^) is strictly positive, namely there 
is K* > such that for any vector v = v{s), s G {1, .., 5*}, 

{v,L^^'^v) = ^L^^\s,s')v{s)v{s') > K*{v,v) (2.6) 

s,s' 

The proof of Theorems 2.2 and 2.3 is given in Appendix C. 



The minimizers satisfy the mean field equation 

p(^)=exp{-/3{^^«-A,}} 



(2.7) 



The Hessian L^'^) has the explicit form: 

L^'\s,s') 



dpsdps 



(k) 



Pp. 



(2.8) 



3 Restricted ensembles 

The purpose of this paper is to study the system in restricted ensembles defined by restricting 
the phase space to particles configurations which are "close to a mean field equilibrium phase" . 
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Unfortunately the requests from the Pirogov-Sinai theory will complicate the picture, but let 
us do it gradually and start by defining notions as local equilibrium and "coarse grained" 
variables, adapted to the present context. 



3.1 Geometrical notions 

We discretize M'^ by introducing cells of size £ > 0, the mesh parameter i will be specified in 
the next paragraph. 

The partition V^^^ 

• P(^), ^ > 0, denotes the partition {cP , x G OJ^} of M'^ into the cubes cf' = {r G M'^ : 
Xi < Ti < .Xj + i?, i = 1, .., d} (vi and xi the cartesian components of r and a;), calling C^^ the 
cube which contains r. 

• A set A is P*^^-* -measurable if it is union of cubes in P^^^ and <^out[-'^] denotes the union of all 
P^^^ cubes in A'^( the complement of A) which are connected to A, two sets being connected if 
their closures have non empty intersection. Analogously, ^f^fA] is the union of all V^'^^ cubes 
in A which are connected to A'^. 

• A function / : M"* ^ M is -measurable if its inverse images are -measurable sets. 

The basic scales 

There are four main lengths in our analysis: Eq <C £-,-y <S ^ ^+,7- More precisely let a+, 
a_ and a verify 

1 . ^ 

->a+>Q;_>a>0 (3.1) 

(the precise meaning of the inequality will become clear in the course of the proofs) , then 

lim ^° „ = lim — . = lim — f^'^ . = 1 (3.2) 

with the additional request that ^+,-y is an integer multiple of 7"^ which is an an integer 
multiple of £-^^ which is an integer multiple of £o- The partition T>^^^ is coarser than V^^'^ 
if each cube of the former is union of cubes of the latter, we will then also say that V^^ ^ is 
finer than V^^\ This happens if and only if ^ is an integer multiple of £', thus is finer 
than T)^^-'i^ which is finer than T>^'^ ^ which is finer than P^^+.t). 
We will need that 

(a^- + a-)d 1 „ ^ 1 . „. 

Eventually we define, for any P^^+'t) -measurable region A, : 

A^A := ^ (3.4) 
where |A| is the volume of the region A, thus Aa is the number of blocks C^^+'^^ inside A. 
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The accuracy parameter C 
Finally, the parameter a in (3.1) is not related to a length, it defines an "accuracy parameter" 

C = 7" (3.5) 

whose role will be specified next. 



3.2 Local equilibrium 

A particles configuration g is a sequence (...r^, Sj....) such that for any compact set A and any 
se{l,..,S}, 

n(x, s) := |g(s) n A| < oo, q{s) = {ri, Si e q : Si = s} (3.6) 
We then associate to any such q the empirical densities 

p(^H^;r,.):=M^^^, se{l,...,S} (3.7) 

as functions on x {l,..,^} and the "local phase indicators" first for any p G L^{W^ x 
{1, .., S}) (p^'^) below as in Theorem 2.3) 



(3.8) 



k if|/ [p(r',5)-pW]|<C,forallsG{l,..,5} 
otherwise 

V. 

and then for any particles configuration q as above, 

r,(C'^)(g;r) = ?7(f'^)(pW(g;-);r) (3.9) 
With ( and l-^-y as in (3.5) and (3.2), we then define 

:= |g : ??('^'^-.^)(g; r) = k, for all r G M*^} (3.10) 

X^''^ is the restricted phase space and the configurations in X^''^ are said to be in local 
equilibrium in the phase k. Their restrictions to a P'^^^'T^-measurable set A is denoted by 
xj^^ and we will study (in the simplest case) the Gibbs measure with Hamiltonian Hx as in 
(1.4) on the phase space restricted to X^'^K As mentioned in the beginning of this section to 
apply Pirogov-Sinai we will need to complicate the picture, by adding a "polymer structure" 
to the phase space and by modifying the Hamiltonian Hx- 



3.3 Polymer configurations 

A polymer is a pair T = (sp(r), Tyr), sp(r), the spatial support of F, is a bounded, connected 
p{^+,7)-measurable region and rjr, its specification, a P^^-'T^-measurable function on sp(r) 
with values in {0, 1, .., 5+ 1}. In the applications of Pirogov-Sinai, T will be contours and rfr 
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not as general as above, to keep it simple wc skip all that sticking to the above definition. We 
tacitly fix in the sequel k E and the corresponding phase space -Y^*^) and define: 



Polymer weights 

The weight of F is a function ^^(r; q), q E X^''\ (dependence on k is not made explicit in w) 
which depends on the restriction of q to 62^^^ [sp{T)] and which satisfies the bound 

sup \w{T; q)\ < e--p°'^'^-,7^r^ = MT)] (g.n) 



Polymer configurations and weights 
We denote by T sequences ...Fj... of polymers with the restriction that any two polymers Fj 
and Tj, i ^ j, arc mutually disconnected (i.e. the closures of their spatial supports do not 
intersect and they are therefore at least at mutual distance i+^-y). The collection of all such 
sequences is denoted by B and B\, A a P^^+'T^-measurable region, the subset of B made by 
sequences whose elements F have all sp(F) in A; B^ subset of Ba with the further request 
that sp(F) is not connected to A'^. If F € i3 is a finite sequence we define its weight as 

wiT-q) = l[wiT;q) (3.12) 

rer 



3.4 The interpolated Hamiltonian 

Pirogov-Sinai applications also require to change the Hamiltonian. Let A be a bounded, 
p(^+,7)-measurable region, q\ G '^^^ , then the "reference Hamiltonian" in A is 

M9A) = E[(E^i'^)-^/3]< E /"H^a;^'^) (3-13) 

s s'^s xeioZ'^nA 

where Xp is the chemical potential introduced in Theorem 2.2, £o is defined in Subsection 3.1, 
in (3.7). 

For any t G [0, 1] we then define the "interpolated Hamiltonian" 

HA,t{qA\qAo) = tHA{qA\qAo) + (1 - t)/iA(gA) (3.14) 

where qA G qA^ G and 

-H"a(9a|9a0 = H{qA U qAc) - H{qAo) (3.15) 

H as in (1.1) with A such that |A — Xp\ < 07^/^. Since -f^A,i(^A|^A<=) = -f^A(<?A|^A':) and 
H A,o{qA\qA'') = ^a(9a), HA,t interpolates between the true and the reference Hamiltonians. 
As we will see in [8], -ffA,t(^A|^A<=) enters in the analysis of the finite volume corrections to 
the pressure, a key step in the implementation of the Pirogov-Sinai strategy. 
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3.5 DLR measures 

The finite volume Gibbs measure in A, A a bounded, P^^+'t) -measurable region, with bound- 
ary condition qj^c , is the following probability on 

dGA(gA,r|gAc) := ^-'^^^ , , du^{qK) (3.16) 

where the free measure dvf^{qjCj is 

I ,^J{l^)d^f^{l^) = y2,~\ ^ I f{ri,si,...,rn,Sn)dri---drn (3.17) 

•^'^A n=0 si,...,s„''^ 

and where the partition function Zf^{qi^c) is the normalization factor which makes the above a 
probability. In (3.16) the boundary conditions only involve particles configurations, to define 
the DLR measures we also need to condition on the outside polymers. 

DLR measures 

Given T G F = (Fi, . . . ), we call F^c the collection of all pairs (sp(Fj) PI A'^, ?7sp(rj)nAc) 
where ??sp{ri)nA'= denotes the restriction of ■qr to sp(F) fl A"^. We then define the probability 
(^G(gA,r|gAc,r^c) on A'f ^ X ^ by 



dGA(gA,r|gAc,f^c) := ^fc^e-''^A,(.A|?Ac) 



n w{T:,q)\dvK{qK) (3.18) 
^rer:sp(r)nA^0 



A probability jx on X^^^ x i3 is DLR if the two properties below hold. 
• it verifies the Peierls bound: for any Li, ..,rjt. 



({r 9 Li} n • • • n {r 9 r^}) < e-^°i^'^-.7(^ri+..+Afr,) (319) 



• for any bounded, P'-^+'^^-mcasurable region A the conditional probability of ^ given that 
the particles configurations in A^ is qt^^ and that F^c = Fyyc is (iGA(^A) FlgA^^) Fyyc) as given 
by (3.18). 

A few remarks on the above definitions: the Gibbs measures dG\{q\, Fj^Ac) satisfy the Peierls 
bound (3.19). Indeed given any Fi, ..,Ffe in such that sp(Fi) is not connected to sp(Fj) 
for any i ^ j, then, for any q\, 

k 

w{T,QA) = {YlH^i,QA)} Yl n ^(r,«A) 

reBA:ri,..rfcer i=i reSA:ri,..rfcer rer,r^ri,j=i,..,fe 

k 

<{[lw{T^,qA)} J2 H^^Qa) 
1=1 reSA 
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and (3.19) follows from (3.11). On the other hand we have not specified all the properties of 
the weights as they arise in the applications (to the continuum Potts model) so that in the 
present context wild things may happen. For instance weights still compatible with (3.11) 
may be such that whenever sp(r) contains (5q^J-^[A], A a bounded, simply connected P^^+.f) 
measurable set, then w{T, q) = unless sp(r) D A. If the weights had such a property then 
there are sequences of finite volume Gibbs measures whose limits are not supported hyVEB. 
Thus a support property like (3.19) is necessary in the present context. 



3.6 Main result 

We fixA;G{l,..,S' + l}, the statements below being valid for any such k and for all 7 small 
enough. We will employ the following notion: {q,T_) agrees with {q' ,]^) in A (A a p(^+>t)- 
measurable set) if all T G F such that the closure of sp(r) intersects A are also in F' and 
viceversa and moreover 

qnA* = q'nA*, A*:=AU{sp(r)U(5S"^[sp(r)]} (3.20) 

rer 



Theorem 3.1. For all 7 small enough there is a unique DLR measure fi and there are con- 
stants ci and C2 such that the following holds. For any bounded, V^^+'^^ -measurable regions A 
and A' D A and any boundary conditions g^c andq'j^ic there is a coupling dQ ofdG\{q\,T\q'j^c) 
and dG\i{q\' ,T\q'J^,c) such that if A is any 'D^^+-^^ -measurable subset of A: 

, . _ dist(A,Ag) 

Q\{ (?A>r') and (gA,,r") agree in A} ) > 1 - cie ^+.7 (3.21) 



3.7 A finite size condition 

The proof of Theorem 3.1 follows the Dobrushin Shlosman approach: we first introduce and 
verify a finite size condition and then prove that this implies uniqueness and exponential de- 
cay. In this subsection we describe the former step. Let A be a D^^+'T^-measurable, connected 
region contained in A* where A* is obtained by taking a cube C € d(^+.t)^ then considering 
A := CuSl'l^ClC] and finally A* = ^U5fui^[A]. All the bounds we wiU write must be uniform 
in such a class. Notice that the diameter of A is > which for 7 small is much larger than 
the interaction range, in this sense A is "large" and we are away from the Dobrushin's high 
temperatures uniqueness scenario. 

Our finite size condition involves only Gibbs measures without polymers: namely the proba- 

(k) - (k) 

bility on Af^ defined for any given q\ G as follows 

g-/3-ffA.t{(?Al9Ac) 

dGliqAlqA.) := -^^j^^-^dMq^) (3.22) 
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We want to compare two such measures with different boundary conditions q'^c and g^c , thus 
introducing the product space x whose elements are denoted by [q'^, q'J^). The finite 
size condition requires that there is a coupling dQ of dG^{qA\q'f^c) and dG^{qA\q'j^c) with the 
property that the event we define below has a "large probability" . 

Notation 

Let m = 2^^ + 2 and 

Cacc — 2c* with c* as in Theorem 5.1 below. Call • — ('eJcC <^iid define 
a partition of M+ into the intervals [0,Cm), [Cm, Cm-i),"-,[C3, C2), [(2,00). 

Definition 3.2. The function Ka{-) and the set 6a( ). 

We denote by 

Ax := Sa;(10~^°^+,-y) n A'', Bx{R) the ball of center x and radius R (3.23) 

Given q'j^c and q'^c, we define the function K\{q'j^c,q'^c',x), x G i-^-^'L^ PI A as follows. 
If Ax = $ then KAiq'/yc ,q'{c;x) = fh+l. 

IfAx^$ and q'j^, D A^ ^ q^, D A^, then KA(g^c, q^c, x) = 0. 

If Ax and q'i,cf^Ax = q'Lf^A^, callb:= max \p^^''''\qAc;r, s) - p'^J'^, then if be 

r<^Ax,s(^{l,..,S} 

[Cm+ij Cm) for some m>2, we set KA{q'\c,q\c', x) = m, otherwise we set KA{Q'Ac,q\c', x) = 0. 

The set Qa{x) = Q A{Q\cq\c', x) , x G l-^yX'^ D A, is defined as the whole space {^A'^aI */ 
K{-;x) = KA{q'j^c,qAc',x) = and otherwise by 

OA{x) = {q'A,qA: q'A n ct '^^ = q'i H ct 

max \p^'-'-'\q'A; x, s) - pW| < Ck(,.)-i| (3-24) 

In section 7.4, we will use Theorem 3.3 below with tt, = 5*^ — 1 and A C A*. Recalling the 
definition of A^a in (3.4), we state: 

Theorem 3.3. For any integer n > there exist 7n > and e„ < 1 such that for all 
7 < 7„ and for any A with Na < n, for any q'j^c and q'J^c as above, there is a coupling 
dQA of dG\{qA\q'A<:) dG^{qA\qAc) such that with K{-;x) = i^A(9A<=' 9a<=' ^) Qa{x) = 
@a{QacQac', x) defined above, 

Qa( n eA(x))>l-e„ (3.25) 



The proof of Theorem 3.3 is given in Part II of this paper. It consists of three parts, in the 
first one we use a step of the renormalization group to describe the marginal of dG^ over 
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the variables {p^^-'''\x,s),x G £-^^Z'^ H A,s G {Ij-.jS'}}. Their distribution is proved to 
be Gibbsian with an effective Hamiltonian at the inverse effective temperature (31- ^-f. In a 
second part we study the ground states of the effective Hamiltonians, proving exponential 
decay from the boundary conditions. In a third and final part we bound the Wasserstein 
distance between the Gibbs measures by approximating the latter to Gaussian distributions 
describing fluctuations around the ground states characterized in the previous step. 



3.8 Disagreement percolation 

The finite size condition established in Theorem 3.3 is used to construct the coupling Q 
of Theorem 3.1. The proof uses the ideas introduced by van der Berg and Maes in their 
disagreement percolation paper, [3]. The proof given in Part III of this paper consists of two 
steps. In the first one we introduce set-valued stopping times, called stopping sets, and prove 
that monotone sequences of stopping sets define couplings of the Gibbs measures and that if 
the sequence stops, then in the last set there is agreement. In the second and last step we 
prove that the probability that the sequence stops late is related to a percolation event which 
is then shown to have exponentially small probability. 



Part II 

The finite size condition 



4 Effective Hamiltonians 

We will use the following notations. 
4.1 General notation for Pcirt II 

• By default in this section A is a connected, P^^+'T^-measurable region contained in A*, see 
Subsection 3.7, and regions in are all P^^-'t) -measurable. To discretize we will use the 
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lattice i-^^l,^. Thus in the sequel l-^^ is the basic mesh. We define 

4\x,y) = \ / J^(r,rO, x,y^a\ 1 = 1.,-, (4.1) 

• The basic variables are the densities pA = {p/\{x-,s) > 0,x G ^-^^Z*^ n A, s G {l,..,^}}, 
A c M'', (by default variables denoted by p are non negative densities). Call the set of 
all PA such that txa := has integer values, so that is the range of values of the 

densities p^^~'''\qA', x, s) when qa G -^a'^ x e i- ^Z*^ n A, s G {1, .., S}; p^^ being defined in 
(3.7). 

• To have lighter notation we will use the label i for a pair (x, s), x G ^-^^Z*^, s G {1, .., S"}, 
writing x(i) = x, s{i) = s if i = {x, s) and sometimes shorthand |i — j| for \x{i) — x{j)\ and 
z G A for x{{) G ^-^-yZ'^ n A. 

• H denotes the Euclidean space of vectors u = {u{i),i G A) with the usual scalar product 
{u, v) = ^2 u{i)v{i)- By an abuse of notation we also denote by H the Hilbert space with A 

i 

above replaced by W^. 



4.2 The effective Hamiltonian 

The effective Hamiltonian H^{pA\qA<'), pA € Xj^\ q\c G is defined by the equality 

e-pei^^HfipAlqAc) ._ f e-^^A.*(«A|9Ac)j,^(^^^) (4 2) 

{p{<ia;-)=pa} 

HA,t as in (3.14), so that is the effective inverse temperature. The Gibbs measure with 

Hamiltonian H]^{p\\q), inverse temperature and free measure the counting measure on 

X^^-* is then the marginal over the variables {p\ G X^'^} of the Gibbs measure dG^{qA\qAc) 
defined in (3.22). 

Since i-^-y = 1-1+°=- and a_ is small, the effective temperature vanishes as 7 — > 0, and the 
analysis of the Gibbs measure becomes intimately related to the study of the ground states 
of H^. This will be the argument of the next section, in this one we determine H'^^. In this 
subsection we describe its main terms and state the main theorem; in the successive ones we 
give the proof. 

The LP term. 

The main contribution to the effective Hamiltonian will be the Lebowitz-Penrose free energy 
functional, the LP term in the title of the paragraph. This is 

i^A(pAlpAc) = t{^ipA,VyPA) + ipA,VyPAc)] - ■^(1a,X(pa)) + (1 - t)(/9WlA, /9a) (4.3) 

where we employ the usual vector notation: if A{i,j) is a matrix, u{i) a vector in H, 

{u, v) = j2 Mi) = J2 ■^')^(-^') (^•^) 

i 3 
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calling 1a the vector 1a(z) = 1 if i G A and = otherwise. In (4.3) 

The normalization is such that is a probability kernel. The term (1a,2'(pa)) in (4.3) is 
"the entropy minus the chemical potential energy" : 

X(pa)(^) = I*iPAii)), T{h) := -6(log h-l) + I3\ph (4.6) 

When t = 1, Fa is just the usual LP free energy and for this reason we call the LP term. 
Notice that if pa(*) = pi(i)lA('')' then the bulk terms of F\ which are proportional to t cancel, 
this will play an important role in the study of the ground states. 

The one body effective potential. 
This term is due to second order terms in the Stirling formula when computing the entropy 
contribution. It has the form: 

K\P^) = ^(lA,log^27r£l,^PA + t[A^ - A]pa) (4.7) 



The many-body effective potential. 

(2) 

This term denoted by Hj^'{p\\q\c)^ takes into account variations of the potential energy 
inside the elementary cells ci^~'^^ G p(^-.7) which have been neglected in the LP term. The 

(2) 

dependence of H)^ on p\ is very simple, it is in fact a polynomial of order < N, N a suitable 
positive integer. The coefficients of the polynomial are described next, they have a simpler 
form once we use Poisson polynomials. We denote by 7rfc(n) = n(n — 1) • • • (n — fc + 1), k G N+, 
n G N+, the Poisson polynomial of order k and, by an abuse of notation we write 

<(p) = rfXn), p=-^ (4.8) 

We shorthand i = (ii, ..,in), n < N, and call n = n{i); i n A / meaning that there is 
ih ^ i such that G A. Given i we denote by k{i) = {k{ii), .., k{in)), with k{ifi) positive 

Hi) 

integers, calling \k{i) \ = ^ A:(z/i). We finally call pho{i) := p''^~'"'\qK'=;i) and denote by p{i) 

h=l 

the function equal to pa(^) and to pA_c{i) when i G A, respectively i G A'^; oq below is a 

(2) 

positive number < 1. Then has the form: 

n(i) 

<^(pA|fAc) =Y. E (7^-,7)"°"®'*aM^),9A^,) n <i^uMih)) (4.9) 

inA^0 fe(i):2<|fe(i)|<Ar h=l 

$ are coefficients which may depend on q\c but only if z n A"^ 7^ 0, in such a case they only 
depend on q\c^i := |^ {qp^c n C'^(7)^^J"" '^^^ main features of the coefficients $ (whose 

dependence on t is not made explicit) is that: 

$(i,fc(i),gAsi) =0 if diam(x(ii),...,x(i„)) > 2Af7-i (4.10) 



16 



and 

X] X] •^*fe^(i)>9Ac,i) < c, for any zo (4.11) 

i3io k(i):2<\k(i)\<N 

where c > is a constant independent of q\c and t. 



Theorem 4.1. For any ao < 1 there are c, N and coefficients $ as above such that for all 7 

small enough 

Hf{pA\qAc) = Fa{pa\pac) + H'f^HpA) + H^aHpaWa^) + RA{pA\qAc) (4.12) 
with the remainder Ra{pa\qA'') = R^^^ + R^^^ 

< cY, i = l,2 (4.13) 
with r = (3 - 5a_ - 2a+)| > (see (4.18) and (4.30)). 



Recall that in this section A is a subset of A* thus |A| < c^+^-y, c a constant, if we wanted larger 
volumes we would have to increase A^, namely to include more body-potentials and longer 
interaction range, the expansion in Theorem 4.1 being highly non uniform in A. The proof 
which follows closely the one in [13] of a similar result, is given in the remaining subsections. 



4.3 Derivation of the LP term 

We fix arbitrarily pA G -^a^\ call nA{i) = i'L^^pA{i), introduce a set of labels C whose elements 
are denoted by ^ = where i = (x,s) € A, £ G {1, ... nA{ii)}', the coordinate functions 

on C are x{^), s(^) and £(^) respectively equal to the first, second and third entry in ^. 
We then define for ^ G A (meaning x{^) G A) the probability measures on A x {1, ..,5'} as 

dpc(r,s) = 1 1 — and call dpA =\ \ dpc, remembering that this measure as well 

as the index set C depend on the initial choice of pA, as this is momentarily fixed we are not 
making it explicit. We obviously have: 



ffdriAii) 



Hfipm = ( rr ^_^) [ e-/^^A,(.Ai?Ac)^p^ (414) 
V4inA(i)!y J 

where on the r.h.s. should be thought of as a ^-labeled configuration of particles (the label 
specifying also the cube where the particle is) which is identified to the integration variable 
relative to the measure dp a- thus the dependence on pA is hidden in the structure of the 
probability dp a- The bracket on the r.h.s. is equal to 

pduAii) 

n^f^ =e'--^'^'^^'^^^ ^(pA)(i) = r;!^(nlog^l,^-logn!), n = nA{i) = £i,^PA{i) 



ieA "^(^) 



(4.15) 
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Then, recalling the Stirling formula: 

n! = n"+^/V"V2^(l + o(^-^^) (4.16) 

we can estimate (1a,5'(pa)) as follows 

(U, S{pa)) = (U, S-pp(pa)) - - m^^^ (4.17) 

where S'^PP(p) = -p(logp - 1) and H^^'^^l is equal to the r.h.s. of (4.7) with t = 0. 
Proof of (4.13) for R^'^\ 

We now show that R^^^ defined in (4.17) above satisfies the bound (4.13): 



ii,^ [5^pp(/3a)(0 - S{pa){{)] = -UAii) (lognA(z) - 1) + lognA(0! 

= ^lognA(i) + logV2^ + f J—^ ] , 

(lA, s^^^ipA)) - (lA, sipA)) = f3H^''"^ + ^ (rf ) 

where we used the fact that nA(i) > c^l,^) since pA G ^a'^- f^om this, we get 

<5iVA(^^yrff' □ (4.18) 

Call H^iqAlQA'') the energy HA^tiQAlQA") defined with J-y replaced by J-y '~' , then i^AC^Al^A^) 
depends only on the densities p^^~'"'\qA',i) and p^^~'^\qA^;i) which in (4.14) are fixed equal 
to pA{i) and pAc(i), hence 

^aC^aI^aO = ^-,'y{t{l{pA,V^PA) + {PA,%PA^) - A(1a,Pa)) + (1 " t)(lA[p^^^ - Xp],PA)} 

(4.19) 

Collecting all the above terms we thus identify in (4.14) 

g-/3i?l,^{Hi"'(pA|5Ac)+R(^)} ^ J ^-P{HAA<lA\qAc)-HA(qA\qAc)}^p^ (4 20) 



4.4 Cluster expansion 

To estimate the r.h.s. of (4.20) we use cluster expansion. Call S a set of unordered pairs (^, ^'), 
^ 7^ then S defines a graph structure {jO.,6) with vertices ( E JC and edges (C^O & £■ We 
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call diagrams the connected sets 9 in (£,<?), 9_ = {Oi, ■■■,0n) their collection. Call and 9dsc 
the spaces of all possible diagrams and of all possible ^ which appear when varying £. Let 

w{e) = / ( n {e-^*{^-'(^(«'^(«'))-^='-v^-'(^(^)'^(^'))> - l})dpf, (4.21) 
then, since dpp^ is a product measure, 

f e-P{HA.ti<lA\qAc)-H^{qA\qAc)}^p^^ Y[w{e) (4.22) 

eeOdsc see 

(4.22) is derived from (4.20) by writing 

^-P{HA,t{iA\qAc)-HAiqA\qAc)} ^ Y[ ^e-(^HViirori')-iZ%Vj{x{i),x{i'))} -1 + 1} 

where the labels ^ include both the particles in A and those of qj\^c outside A. After expanding 
the product we then get (4.22), details are omitted. 

The basic condition for cluster expansion which we have in the present context, involves the 
elementary diagrams namely 9 = (^, (') and states that given any a > 

\w[{^, ^')) |7~°-+" < 1, for any 7 small enough (4.23) 

(4.23) is proved by observing that the densities pA{i) are bounded and that (4.21) yields for 

l»((«.0)l<^^7^-.)l,,,,,<-.,,,.,-.,,,,_. (4.24) 

"Cluster expansion" then applies for any 7 small enough and the following holds (for any 

'-A 



PA e 4'^ 



Notation. We give O a graph structure by calling vertices the diagrams 9 & Q and edges the 
pairs 9 and 9' which have non empty intersection, as sets in C 

Denote by m(0), G 0, positive, integer valued functions, calling m{9) "the multiplicity" of 
9. We restrict to m & M where 

m G if and only if sp(m) := : ^ G 9,m{9) > 0} is a connected set (4-25) 

and shorthand ^ G m when ^ G sp(m). 

Cluster expansion tells us that given any oq < 1 for all 7 small enough there are coefficients 

Lo{m), m & Ai, such that 

logZ{{w{-)}) :=log{ llw{9)} = '^M (4-26) 



and, for any ^ & C, 



\uj{m)\{ W (7£_^)«ol^l''<ig"^(^)l} < 1 (4.27) 

meM:m3£, 9:m{e)>0 
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where |^|edg is the number of edges in 9. The coefficients u{m) have the following explicit 
expression: 

uj{m) = Cm n (^-28) 

e:m{e)>o 

where thinking of Z({tt;(-)}) in (4.26) as a function of the weights {w{9), 9 G O}, 



n rn(9)\ ^ n 



Qm{e) 



0:mi0)>O ^ ' e:m(6)>Q^^^^> 



} log^A(w;(-)) 



«)(6l)=0 



(4.29) 



(Cfn being bounded coefficients independent of A). As said, all the above follows from the 
general theory (of cluster expansion) using the condition (4.23), see for instance [16]. 



4.5 Identification of the many body potential 

We will next use (4.27) to truncate the sum in (4.26) identifying the remainder with the term 
and recognizing in the finite sum the Hamiltonian i?^^^(y9A|9A'=); for this we will use the 
explicit representation of the terms of the expansion provided by (4.28)-(4.29). 
Calling \m\ = ^ |6'|edgm(6'), by (4.27), for any N > 0, 
eee 

meM:\m\>N ^eC m(^M:m3(,\m\>N 

< |£|(7^-,^)»°^ ^ |a;(m)|{ J] (t^-.t))""""'^^^"^^'^'} < I-^Kt^-.t)""^ 

meM:m3i 6»:m(0)>O 

Since A C A*, there is c > such that \jC\ < c£+,-y and we can then choose N so large that 

-/3£li?(2) := i^(rn), | ^ uj{m)\ < £1"^^^ (4.30) 

meM:\m\>N meM:\m\>N 

thus (4.13) is satisfied and 

meM:\m\<N 

The dependence on pA is hidden in the space 9, on which the functions m are defined. 
Theorem 4.1 will be proved once we show that the r.h.s. of (4.31) can be written as the r.h.s. 
of (4.9). 

We rewrite the r.h.s. of (4.31) by first summing over all m in "the same equivalence class" 
and then summing over all equivalence classes. Before defining the equivalence m ~ m' we 

observe that if is a one to one map of C onto itself, then ^ extends naturally to a map of Q 
onto itself by letting tp{9) be the diagram with vertices V'(O) ^ ^ ^) and edges (V'(C)) V'(^'))) 
(^, ^') the edges of 9. We then call m ^ m' if there is a one to one map from C onto C such 
that • x((/)(0) = x{0, s(<^(6) = s(0 for all • m'{(f){9)) = m{9) for all G 6. 



20 



Calling [m] the equivalence class of m, i.e. the set of all m' : m' ~ m, we define the average 
weight 



card(H) ^^^j 

Notice that if sp(m) consists only of ^ such that x(^) G A then a;(m) = u{m') = oo*{m) for 
all m! G [m]. If instead there are labels ^ in sp(m) such that x(^) G A"^ then a;*(m) is a non 
trivial average. Actually the averages involve the labels (. in each triple (x, s,^), x G A^, with 
m{x, s,i) > 0. Calling K{i; m) the number of ^ G m such that i(^) = i, 

c&rd{[m]) = Yl7TK{i;m)inii}) (4-33) 

i 

where 7rfc(n) is the Poisson polynomial and n{i) = p{i)£'L ,y. We then have 

-/3£iH'-'\pA\qA^) := J2 ^*{m){ll7TKi^■,m){nm (4.34) 
[m],\m\<N i 

We next interchange the sums: for any sequence K{i) G N+, Yk{i) < N, let 

i 

^{K{.)):=£Z', u;*{m)l[et'^f (4.35) 

[m],m:K(-;m)=K(-) i 

then 

-PH('\pA\qAc) := J2 ^{KmU^'-^^^'^m^inm (4-36) 

thus identifying $ in Theorem 4.1 in terms of 

*(K(.)) = (7^-,^)'^«l^®lf(^,i^(^),gA^,) (4.37) 

recalling the remark before (4.33), indeed the l.h.s. depends on q\<: only via gA=,i- 

Of course we still need to prove that the function $ defined via (4.37) satisfies the bounds 

stated in (4.10)-(4.11). Since the coefficients Cm in (4.28), are bounded, say 

max \Cm\ < Cat (4.38) 

m:\m\<N 

we just need to bound |w(^)|. The definition of w{9) involves product of terms 'u^((^,^')) foi" 
each edge of the diagram which wc bound using (4.24). The bound obtained in this way is 
the same for all m' G [m] so that the bound for uj*(m) is the same as for uj(m). To fix up the 
combinatorics, we proceed as follows. For any m we define a graph structure G{m) on sp(m) 
introducing a node for each element ^ of sp(m) which is then given the label i = {x{$,), s(^)), 
thus different nodes may have the same label. Edges in G{m) are the union of all the edges 
present in all the diagrams 6 such that m{6) > 0. Each edge is then given a multiplicity 
equal to the sum of all m{9) over the diagrams 9 which contain the given edge. With this 
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definition any m! G [m] gives rise to the same G{m) as we are only recording the coordinates 

x(^) and ,s(^) of ^. 

To proceed witli tlie bound we assign a "weight" £l ^ to any node in G{m). Having (4.24)in 

mind, we assign to each edge a weight ( cy'^i'yi- 'y)l le .) ie ^) 1 , where p the 
^ & \^ ? V/ ,7; dist(c;^-'',r/j-,'^')<7-i; 

multiphcity of the edge. We have thus assigned a weight W{G{m)) to G{m) equal to the 

product of the weights of its nodes and of its edges and, with reference to (4.35) and recalling 

(4.38) 

\^{K{-))\<CNrJ.^ W{G{m)) (4.39) 

[m] ,m:K{-;ra)=K{-) 

Recalling that K{i] m) is the number of G m such that i(^) = i, K{i; m) is also the number 
of nodes in G{m) with label i. Thus, calling K{i, G) the number of nodes in G with label i, 
i = {i,ie G}, and K{i, G) = {K{i, G),i G i}, 

\^{K{i))\<CNr'^^ W{G) (4.40) 

G:K{i;G)=K{i} 

(4.37) then yields 

Mi, K{i),qA^,i)\ < CArr;!^(7^-,7)-"°l'^®l J2 ^(^) (^-^l) 

G:K(i;G)=K(i) 

By (4.35) the terms to consider have i such that ^^ivr(z) < A^. Then ^{i, K{i),q\c^i) = if 

diam(x) > 2j~^N, x being the sites appearing in i, because the weight of the edges in G are 
proportional to 1 i- 

To prove (4.11) we fix io and restrict the sum in (4.41) to G : K{iQ] G) > 0. For each such G 
we can then define a tree structure rjg(m) in G{m) with root zq, a first generation made by 
all nodes connected to the root, second generation made by the nodes connected to those of 
the first generation and so forth. To recover the original graph we may also have to add edges 
connecting individuals of the same generation and also attribute to each edge its multiplicity, 
as explained earlier. We then have 

l.h.s. of (4.11) < Y E ^:i(7^-,7)~"°'''®' E ^(^^o) (4-42) 

iBio K{i):\K{{)\<N Ti^:K{v,Ti^)=K{i) 

Define a new weight W*{T) by changing the weights of the edges into 

^c(7l_,-y)^~"°7'^l|j._3;/|<2-y-i j , p the multiplicity of the edge 
while the weights of the node are unchanged. Then 

l.h.s. of (4.11) < rf,^ E ^*(^*o) (4.43) 

i3io K{i):\K{i)\<N Ti^ ■.Kii;TiJ=K{{) 
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The weight of the root of the tree cancels with the prefactor ^I^- We upper bound the sum 
on the r.h.s. if wc regard a multiple edge with multiplicity k as k distinct edges originating 
from a same node and also regard edges between nodes in the same generation as edges into 
the next generation (thus dropping the constraint that the arrival node is the same as the 
arrival node of another edge), each node added in this way getting an extra weight In 
this way we have an independent branching and since 

lim^(7£_,^)"07'^l|,,|<2^-i£l = 

x' 

we then get (4.11), details are omitted. Theorem 4.1 is proved. □ 



5 Ground states of the effective Hamiltonian 

In this section we study the ground states of the main term in the effective Hamiltonian 
H]^{p\\q\c), which, with reference to (4.12), is 

/(pa; QA'^) ■= Hf{pA\qAc) - -RA(pAlgAc) (5.1) 

While originally pA = (^pA{i),i = {x,s),x G i-^^Z'^ H A,s G {l,..,^}) G defined in 

Subsection 4.1, it is convenient here to extend the range of values of pA{i) to an interval of 
the real line. We thus call 

) = [pa : Pa{x, s) e [p« - C, pi') + C], G £_,^Z'^ n A, V5 G {1, .., S}] 

The ground states in the title are then the minimizers of /(pa;9a=) as a function on Yj^^ 
with gA<= regarded as a parameter. 



Let Ka{x) = KA{q'Ac, Qac', x) be the function defined as Ka{x) in Definition 3.2 but with the 
set Ax in (3.23) replaced with the set 

= 5^(10-3%,^) nA^ (5.2) 

Our main result is the following theorem: 



Theorem 5.1. There are c* and Co positive such that for any ao < 1 and for all^ small enough 
the following holds. For any qa^ G xj^c there is a unique minimizer pA o/ {/(pa; ^A=), PA G 
yj^^}. Let K{x) X G £_^-yZ'^n A, be as above and p'j^ and p'^ the minimizers with q'j^c and q'^c, 
then for any s G {1, .., S}: 

• (i) Ifk{x) > 0, \p'a{x,s) - pI{x,s)\ < ce-iO""'(T^+.7)<^. 
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(ii) Ifk{x) = m > 0, s) - pf \ < c*(Cm + (7^- ,7)"° + 6"^° '°(t^+.7)'^), with same 

bound for p'ji^{x, s). 



Existence of a minimizer follows from / being a smooth function on a compact set of the 
Euclidean space. Uniqueness and exponential decay are more difficult and the proof will take 
the whole section. The basic ingredient is that D^/ (the Hessian matrix of the derivatives 
w.r.t. the variables PA{i)) computed on the minimizer in the constraint space Y^'^^ is positive 
and "quasi diagonal" , which would then give the required uniqueness and exponential decay 
if we had Df = 0. This is however not necessarily the case because the minimum could be 
reached on the boundaries of the domain of definition, which, on the other hand, is necessary 
to ensure convexity. We will solve the problem by relaxing the constraint and then studying 
the limit when the cutoff is reconstructed. 



5.1 Extra notation and definitions 

The basic notation are those established in Subsection 4.1, here we add a few new ones specific 
to this section: 

• We will write /{pa; qac) = F{pa; pa'=) + g{pA; qA<=) where, recalling (4.12), 

g{pA; QA^) = H^a\pa) + i^f (paI^aO (5-3) 

• To evidentiate some of the variables in pA, say those in A C A, we write pA = {PA, Pa\a)j 
where pA and Pa\a the restrictions of pA to A and respectively to A \ A. 

(k) (k) (k) 

• It will be convenient to relax the constraint pA £ 5^ by enlarging into VFy^ 

WP = {pa ■ PAix, s) G [pi''^ - b, pf^ + 6], Vx G i-^^l^ n A, Vs e {1, .., 5}} (5.4) 

Wpiki) ^ p{k2)\\ 

where b := min — has been chosen such that 

k\^k2 2 

We then introduce a cutoff parameter e G (0, 1) (which will eventually vanish), call (a)+ = 
ala>o, (o)- = ala<o and define for any e > 0, the function /g on as 

feiPA-, qac) := /(pa; qac) + ^ E ({(/^a(^) - iP% + 

ieA 



24 



• Since / [/g] is a continuous function of pA which varies on a compact set, it has a minimizer 
denoted by p\ [pA,^], and we will later see that this minimizer is unique. We call p its extension 
to the whole l^^^U^ x {1, .., S*}, by setting p = pA^ on A'^. Here pA= is the density associated 
to gA= via (3.7) with ^ = thus p of course depends on q^c 

• For any P^^-'T^-measurable set B we write for any differentiable and P^^-'T^-measurable 
function '^(p) 

Db^ = { xii) e n b] (5.6) 



5.2 A- priori estimates 

In this subsection we prove some a-priori bounds on PA,e{i)- When e > we loose the bound 

Ipa(^) — — C valid at e = but, as we will see, we have the great simplification to 

know that for e small enough, minimizers are critical points, thus satisfying D\f^ = 0, and 

\PAAi)-p%\<^C- 



Lemma 5.2. There is a constant c > such that for all e > and for any minimizer 
PA,e £ wj^^ of fe the following holds: for all x G l-^^TL'^ fl A and all s G {1, .., S}, 

|pA,e(x,.)-pi'=)|<C + c(^)'^/V/^ (5.7) 

(k) 

In particular, if C < b/2 then for all e > small enough, any minimizer pA,e G W^a "/ 



also a critical point. 
Proof. We denote by 



V'(pa) = E{(PA,e(i) - [p% + c])+}' + mAi) - [p% - c])-r 



ieA 



Then for aU PA £Wl^\ 



^V'(pA,e) < /(pa; QA'^) - /(pA,e; QA-) + ^V'(Pa) 



(k) 

and since ^ vanishes on : 



;J-V'(pA,e) < inf /(pa; qac) - /(PA,6; ^a^) 



and, calling cf)' = miri /(pa; qa^), = min /(pa; ^a-^) 

1 



^^^{pA,e) <<P -(p' 
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and in conclusion 

\PaAx,s) - pW| < [H<f>' - ^"))'^' + C (5.8) 
and (5.7) follows because (f)' and 0" are bounded proportionally to the cardinality of {x : x G 

£_,^Z'^nA}. 

By choosing e so small that ( + c7"*^°++""^'^/^e^/^ < 2^ < 6, we conclude that p\^e is in the 

(k) 

interior of VF^ and is thus a critical point. □ 



Lemma 5.3. p\^e converges by subsequences and any limit point p\ is a minimizer of f. 

Proof. Convergence by subsequences follows from compactness and by (5.7) any limit point 
PA is in yI''\ Now for any pA € Y^''\ we get /(^a) = /^(pa) > fe{pA,e) > /(PA,e) and by 
taking e — > along a convergent subsequence /(pa) > /(pa)- D 



A minimizer pA of / is not necessarily a critical point, i.e. D\f = 0, the equality may fail if 
the minimizer is on the boundary of the constraint. In such a case however, the gradient if 
different from zero "must be directed along the normal pointing toward the interior" . 

Lemma 5.4. Any minimizer p\ of {f{pA,QA<'), pA € ^a^^} *^ critical point" in the 

following sense: 



• If for some i G A, \pA{i) — Ps(i)\ < C (strictly!), then 



d 



If instead pA{i) = P^'(\ i then 



dpA{i, 

(i) 



-/(pA,gAc) = (5.9) 



d 

/(pA,(ZA'=) < 0, respectively > (5.10) 



dpA{i) 



5.3 Convexity and uniqueness 

Convexity is a key ingredient in our analysis: 

Theorem 5.5. Given any k G (0, k*) (k* as in (2.6)), for all 7 small enough the following 
holds. Let pA G wj;^^ be such that \pA{i) — P^^^i^l ^ then the matrix A := D\f^{pA, ^a<=) is 
strictly positive, as an operator on Ti, namely (recall the definitions in Subsection 4-1) 

(^u, Au) > K.{u,u), forallu^TL (5.11) 



26 



Same inequality holds when e = 0. 

Proof. Recalling (5.3) and denoting by pj^ below the diagonal matrix with entries pK{i)~^ 

{u, Au) = t{u, V^u) + ^{u, pl^u) + («, [Dlg\u) + {u, [Dl{f, - f)]u) 

and get a lower bound by dropping the last term thus reducing the proof to the case e = 0. 
Extend u and A as equal to outside A and set 

U{x,s)=£i^^ 4-'''\x,y)u{y,s), x e 

where J^f^ is defined in (4.1). Then, 

{u,Au)>tJ2 Yl Uix,s)Uix,s') + ^{u,pl\) + {u,[Dlg]u) 

- [p^^^-^*P{x,sf + ^{u,p]^\) + iu,[Dlg]u) 

recalling(2.8), by (2.6) the curly bracket is non negative as well as — tt^ — k*- 
Since for each s 

U{x,sf< «(^'«)' 

then 
Thus 

Recalling (A.2), (A.3) and using (4.9)-(4.11) we get 
Thus 

(n, [Dlg]u) < [^i.,^f''^{u,u) 
(5.11) is then proved recalling the assumption |pA(i) — p[''(j^\ < 4^. 



□ 
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Theorem 5.6. Given any k G (0,/t*) (k* as in (2.6)), for all 7 small enough the following 
holds. Let /5A,e be a minimizer of /g and for e = of f, then for both e > small enough and 
e = 

fe{pA, > fe{pA,e, QA-) + ^ {PA - PA,e, PA " PA,e) (5.12) 

for all PA such that \pA(i) — P^%\ < 2^ for all i G A. (5.12) remains valid if pAe is a critical 
point, D\f^ = 0, and \pA,e — P^''^ < 2( as well as when e = and p\fi a "critical point" of 
f in the sense of Lemma 5.4- 



Proof. We interpolate by setting pa{0) = Op\ + (1 — d)pA,e, d £ [0, 1], then calling V'e(6*) '■= 
fe{pA{0)AA<^) we have 

V'e(l) - ^-6(0) = C {DM0),pA-pA,e) 

Jo 

= [ [ {Dli^,ie'){pA-pA,e},PA-pA,e) + {DAM^),PA-pA,e) 

Jo Jo 

By (5.7) for e > small enough and for e = as well, \pA{d) — P^''^\ < so that by (5.11) 

Jo Jo ^^'^^^'^^"'^^^ ~ PA,e},PA - pA,e) > ^ {pA - PA,e, PA - PA,e) 

Moreover (Da^£(0),pa — pA,e) > 0- In fact, if e > and pA,e is a minimizer of /g, by Lemma 
5.2 (for e > small enough) pA^e is also a critical point and Z?AV'e(0) =0. If e = and pA 
a minimizer of / then by Lemma 5.4, (-DaV'o(O), pA — Pa) > which, for the same reason, 
holds if PA is a critical point of / in the sense of Lemma 5.4. 

□ 



Corollary 5.7. For any 7 and e > small enough the minimizer of f^ is unique, same holds 
at e = for f. For e > (and small enough) there is a unique critical point in the space 
{|PA — P^*^^! < 2^}; such a critical point minimizes /g. Analogously, when e = there is 
a unique critical point in the sense of Lemma 5.4- Such a critical point minimizes f. The 
minimizer of f^, e > 0, converges as e ^ to the minimizer of f. 



Proof. From Lemma 5.2 it follows that any minimizer pA.c of is also a critical point and 
verifies (5.7), so for all x G l^^^ll^ n A and all s G {1,..,5}, |pA,e(a;,s) — p^s \ < 2^ and 
we can apply Theorem 5.6 to the matrix D\fe{pA,e'-, QA'=) ■ If we assume that there are two 
minimizers, then (5.12) gives a contradiction. The proofs in the case e = follows by using 
Lemma 5.3. □ 
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5.4 Perfect boundeiry conditions 

In this subsection we restrict to "perfect boundary conditions", by this meaning that we 
study 

FHpa; Qa^) = FA{pA\p^''hAc) + g{pA; qa^) (5.13) 

namely wc replace in the LP term of the effective Hamiltonian, see (4.12), by the mean 
field equilibrium value, ff^ is then defined by adding to /P^ the term fe — f given by (5.5). 
All the previous considerations obviously apply to /p^ and fi^. 



Theorem 5.8. For any 7 small enough and for all e > small enough, the minimizer 
of ff^ minimizes f^^ as well and it is such that 

IPaA'^ - P%\ < c(7^-,t)"°, for alii e A (5.14) 

c > a constant. 



Proof. Since p^^ is a minimizer of ff^, DAfF^ipA^ = 0. Then if (5.14) holds, DAf^^{pAe) ~ 

DAft^{PAe) ~ ^ ^^'^ Corollary 5.7 fP^^ is a minimizer of /p^ We thus have only to prove 
(5.14) for all e > small enough. Consider first the simplified problem with g = Q. 

Case g = 

Recalling (4.3), if -Da-Fa(pa|p*^*^^1ac) = 0, then by an explicit computation, for all i G A, 

PA(^) = exp{-/3[ J2 iV'yiiJ)pij) + i^-t)P%-h\} (5-15) 

where p(j) = pa(j) if i G A and = p^^j-^ if j G A'^. pA{i) = is a solution of (5.15) and 

therefore also a solution of DAff^ = (with g = 0). By Corollary 5.7 it is then the unique 
minimizer of ff^ and (5.14) is proved (for g = {)). 

Proof of (5.14). 

Call 

fe,e{p^) = ^a(pa|p^'=)1aO + ^5(pa; ^a^,*) + {fe - f) 

9 G [0, 1] ; for all e > small enough denote by pA,efi the minimizer of f^fi, so that DAfe,e{PA,e,e) 
0. Suppose that 

^^j^'^ exists for all 6 G [0, 1] and depends continuously on 9 (5.16) 
du 

Obviously PA,e,i = PA,e while pA,e,o = P^'^^^A because of the above analysis with g = 0. Then 

PA,e,e = P^lA+ (5.17) 



d9 
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On the other hand by differentiating I?A/e,e(pA^g^e) = we get 

DUAPA,e,e)^^ = -DAg{pA,e,e) (5.18) 

By Lemma 5.2 and Theorem 5.5 for all e > small enough, -D^/e,e(/OA,e,e) is symmetric 
and positive definite, then by Theorem A.3 the inverse {D\f^fi{p\fi))~^ is well defined and 
bounded as an operator on and we thus get from (5.18) 

l^^l < c\\DAg{pK,e,e)\\oo < c'{li-,,r (5.19) 

which by (5.17) yields (5.14). (5.19) also implies that \pA^e,e - P^^'^IaI < c'(7^_,-y)"o. Notice 
that (5.19) implies (5.16), but unfortunately the argument is circular as it started by suppos- 
ing the validity of (5.16). To avoid the impasse we start from the equation in the unknown 

Dlf,,0{pA)uA = -DAg{pA) (5.20) 

where pA is considered as a "known term" such that \pA{i) — pI'(1^\ < '^C for all z G A. From 
what said before, (5.20) has a unique solution called Pa(^|pa) and 

\pA{i\pA)\ < c(7£_,-y)''°, for alH G A (5.21) 

Since pa{-\pa) is Lipschitz in pA (we omit the details) the ordinary differential equation 

= Pa(>a(^)), PA(0)=p(^)lA (5.22) 



de 

has a unique solution pa{0)- Then, by (5.20), 

^DAfeApAi^)) = 0> and hence L'A/.,e(PA(-; 0)) = DAfeM'^'^W = (5.23) 

Since \pA{e) - p^^^^aI < d{-^^-^^Y\ DAh,e{pA{-,0)) = as well, hence by Corollary 5.7, 
Pa{-',G) = PA,e,e{') and by (5.22) it is differentiable with continuous derivative. (5.16) thus 
holds and the theorem proved. □ 



5.5 Exponential decay 

This subsection concludes our analysis with the following main theorem, Theorem 5.1 will 
be proved in the Subsection 5.6 as a corollary, taking Af as a neighborhood of x in A'^ and 
A| = A^\Af. 
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Theorem 5.9. There are Co and c positive such that the following holds. Let p'^ and p"f^ he the 
minimizers of f{pA, q'\c), respectively f{pA, <i'\a), with q\c,q\c € -^Ac Then for any partition 
of into two T>^^-'^^ -measurable sets and A^, 

IpaW - Mi)] < c(mm{l,»^^,. ^;max ((7^-,^)"° + |/--)(^y(c; j) - P^'-'^HqU 

^ 1 1 jfc^i 

+ ^ V^GA (5.24) 

Proof. We follow the interpolation strategy used in the proof of Theorem 5.8. To this end we 
separate the "interaction part" in writing /g = /° + where = f^{pA) is independent 
of the boundary conditions while 

f}{pA, QA'^) = t{pA, VyPAo) + gi{pA, qaO (5.25) 

where gi is given by the r.h.s of (4.9) with the sum over i restricted to the set i H A'^ 7^ 0. 
We then interpolate between the two boundary conditions 

feAP^) ■■= feiPA) + dfeiPA, 9ac) + (1 - e)f}ipA, q'ac), € [0, 1] (5.26) 

The analysis done in the previous subsections, applies to fe,e{P^) ^ Thus the minimizer 

PA.efi of fe,t is unique, is a critical point, namely DAfe,e{PA.e,e) = ^^'^ satisfies for all 

X £ i-^^Z'^nA and ah s G {l,..,^}, \pA,eA^^^) " Ps'^^l ^ 2C. 

We can apply the same proof as the one given in Theorem 5.8. In fact by Theorem 5.5, for all 
e > small enough, and for all pA such that |pA(a;, s) — pi''^\ < 2^, we have that Dj^fgApA) 
is symmetric and positive definite, then by Theorem A.3 the inverse (Dj^feApA))"^ is well 
defined and bounded as an operator on L°°. Thus the equation 

AdS^ / dDAfe,e{P^) iKr,^\ 

{DAfeAPA))uA = (5.27) 

has a unique solution that we call ua{-,Pa) that is Lipschitz in pA- This implies that the 
equation 

dpA,e,e I ^ 

= UA{-,PA,e,e), PA,e,0 = PA,e,0 

has a unique solution that coincides with the minimizer PA,e,e- Thus PA.efi is differentiable in 
and dpA.e,e/d9 satisfies 

m2f /A xx^^PA^_ dDAfe,e{PK^,e) i^^c.-. 

(^A/(?,e(PA,e,(?))— ^ (5.28) 

By Corollary 5.7, PA.efi converges by subsequences as e — > to a limit pAfi which minimizes 
fe, so that 

|pX(^)-/5U^)l<lim^'|^^^| (5.29) 

We now estimate uniformly in e and 6 to prove (5.24) as a consequence of (5.29). 

dv 
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we let 



Equations for dpA.e,d/d0- 
" ^ ^ -^^A/e',e(PA.e,e) ^ := £>A/e,e(PA.e,e) (5.30) 



so that (5.28) becomes 
We also define: 



Au = V 



Ao:= Dlfefl{pA.cfi), a:=A-Ao (5.31) 
a is a diagonal matrix whose diagonal elements are 

a{i) := 36-1 ({(^^^ ,(^) _ [pW + ^])^|2 + {(^^^ ,(^) _ [pW _ ^])_}2^ (5 32) 

To distinguish among large and non large (called small) values of a{i), we introduce a large 
positive number b which will be specified later and, calling H the Hilbert space of vectors 
u = (u{i),i G A), 

G = {(?;) : a(i) >b}, nG = {u€n: u{i) = 0, for all i G G"} (5.33) 

Let Q be the orthogonal projection on Hg and P = 1 — Q, thus Q selects the sites where a 
is large and P those where it is small. 

Our strategy will be the following: rewrite Pu, Qu as linear expressions of Pv, Qv to get 
bounds on Pu, Qu (and therefore on u) using knowledge on v. 

Rewriting Pu,Qu in terms of Pv,Qv. 
Since the matrices a, P, Q are diagonal they commute, giving for instance QaP = aPQ = 0, 
i.e.: 

QAP = QAqP, (5.34) 

and symetrically: 

PAQ = PAoQ. (5.35) 
Using Q"^ = Q together with (5.34) we get: 

QAQQu = QAQu = QA{u - Pu) 
QAQQu = Qv- QAPu 

Qu = {QAQy^iQv - QAoPu} (5.36) 

where QAQ is invertible on the range of Q since yl is a positive matrix. 
Using P^ = P together with (5.36) and (5.35)we get: 

PAPu + PAQu = PAu = Pv 
PAPu + PAo{QAQy^{Qv - QAqPu} = Pv 

(^PAP - PAoiQAQy^QAojPu = Pv - PAo{QAQ)-^Qv 
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Let 

B = PAP - PAo{QAQ)-^QAo (5.37) 
so that if B is invertible on the range of P (as we will prove) , then 

Pu = B-^{Pv - PAo{QAQ)-'^Qv} (5.38) 



A decomposition of v. 
Recalling (5.30) and (5.26), after expanding the Poisson polynomials in (4.9) we get, 

v{i) = -tY,V,{iJ){pUj)~pUj)) 

-dn{ii, ki^,.., in, kir,]q'K-]t)p'{iif'i-^ ■ ■ ■ p'iinf'") (5.39) 

where p"{i) = p'{i) = PA,e,e{'>') if a:(i) G A and p"{i) = p'l^c{i), p'{i) = p'/^c{i) when x(i) G A'^. 
The coefficients dn satisfy the same bounds as the coefficients $ of (4.9) (with maybe a 
different constant). 

Shorthand by {xj} the sites in {x{ii), ..,x{in)} which are in A*^, noticing that by definition 

of gi there are not terms with {xj} = 0. 

We then call v^^^ the sum of —t ^ V-y{i,j) (j^Ac{j) — p'\c{j)^ minus the second sum on the 

r.h.s. of (5.39) restricted to sets (n,...,z„) such that: {xj} ^ and any Xj G {xj} is either 
in Kl or q'' = g' _ , = C^'^\ (or both), v^"^^ := v - v^^\ 

By linearity u = u^^^ + u^^^ where u^^^ and u^^^ are defined with v replaced by v^^^ and v^^^ 
and we will bound differently u^^^ and u^^^ using || • ||oo norms for the former and || • || norms 
for the latter. 

Bounds on u^^^ . 
By Theorem A.l if 6 is large enough and c > ||^o||) 

\\PAo{QAQ)-^QAo\\ <^=:6, \\PAo{QAQ)-^QAo\\^ < ^e^^^' (5.40) 
Moreover by (A. 5) 

supV |5(z, j)|e^l^-^l < sup y \A{i,j)\e^\'-^\ + < c"'b =: a (5.41) 

Then applying Theorem A.2,A.3 with B as in (5.37) and Ri = PAo{QAQ)-^QAo, B is 
invertible and there is a constant c > such that ||S~^||oo < c. Therefore there is a new 
constant c such that 

\Pu^^\{)\ < cmax \v^^\j)\ (5.42) 
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If q'l^c = ^c, v^^^ = and u^^^ = as well, let us then suppose q'J^c 7^ g^c. Then (5.39) yields 

\Pu('\i)\ < c(m£C| IpUj) - PAc(j)l + (7^-,7)"°) (5.43) 

To bound \Qu^^\i)\ we go back to (5.36), the same arguments used before prove that 
||((3^<3)~^||oo < c as well, so that \Qu^^\i)\ is bounded as on the r.h.s. of (5.43) (with a 
new constant c) and |u(^)(i)| is therefore bounded as the first term on the r.h.s. (5.24), we 
will prove next that |u(^^(i)| is bounded as the second term on the r.h.s. (5.24) which will 
then be proved. 

Bounds on u^^^ . 

Recalling the definition of v^'^^ 

\v'^'Ki)\< ^i^^(i,j)V (5.44) 
where ""^^ K^{i, j) < ck and K^{i,j) = if \x{i) — x{j)\ > c'7~^, c and c' suitable constants. 

i 

By Theorem A. 2 

\B-\i,j)\ < (- + ^)exp| - k' = k-6, S as in (5A0) (5.45) 

a K L a + k' i 

By (5.45) and (5.44), calling c" = 1/a + and u = k' /{a + k'), 

|S-ip^;(2)(i)| < _{ci^c"e"''^}e-'^^l^«-^(^')l (5.46) 



CT 'C~ 

^2 



By (A.5) 

Y^\{QAQ)-\z,j)\e^\^-^\<^ (5.47) 

i 

and since Ao{i,j) = if |z - j\ > c'^~^ and ^ \Ao{i,j)\ < caq, 

i 

\B-'PAo{QAQ)-'Qv(^\i)\<Y,Yl E {c"e-Tl-«--(^')lcAoe'^''^} 
xe-7N(/')-^(/)l|(Q^g)-i(/,/)|eTl^(/')-^{/)lif^(/^/'') 
< {^'cA.e'^^} E h"- _ e-^l-«-(^"")le-' (^) ck 

Thus supposing w < 1, we get from (5.38) 

\Pu^'^\i)\ < Y h" ^P' ce-'^Tl^W-^^^)! (5.48) 
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To bound Qu^^) (recall (5.36)) we use (5.47) to get 

\{QAQ)-^Qv^^){i)\ < Y: ^^e-^l-W-WI (5.49) 



c. c. 

^2 



while, using (5.48) and (5.47), 



hence 



|Qu(2)(i)| < - ^ e-^l-»--(^)l (5.50) 

□ 



5.6 Proof of Theorem 5.1 

The proof is a corollary of Theorem 5.9. Indeed given any x € l^^^U^ H A, call the union 

of all &y"-''\ y G i^^yZ'^ n (A^ n B^{10~^^i+^y)) . Then if k{x) > 0, same notation as in 
Theorem 5.9, q'^c = q'j^c and by (5.24) we are reduced to a sum over j G A2. We split the 
exponent —^ijj\x{i) — x{j)\ into two equal terms and get 

\pI{x,s)-p'j,{x,s)\ < c{e-(-/2)7[«+,.-^-,.]}{^ e-(-/2)7k--0)l| 

< ^/g-(u;/2)7[10-30^+,^-^_,,] ^5 5-^^ 

The exponent ci) in Theorem 5.1 is thus going to be half the uj of Theorem 5.9. Using Theorem 
5.9 with p'j^c replaced by p^'^^Iac, and calling p^ the corresponding minimizer, 

\g{x, s) - Mx, s)\ < c'e-(-/2)Mio-^°^+.-^-,.] + (ci(7^-,,)"° + Cm) 
and using (5.14) 

\p'a{x, s) - pW| < c'e-(-/2)^[io-"^+--^--l + ([ci + c](7^_,,)«o + C„.) 
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6 Local Couplings 



In this section wc prove Theorem 3.3, thus we fix a region A, union of a finite number N\, of 
cubes of and two boundary conditions qi^\c G xj^J , i = 1,2. Wc also fix a t G (0, 1] and 

we consider the two Gibbs measures dG^{qA\Qi,A'') « = 1,2 defined in (3.22) and with state 
space The aim is to construct a coupling Q\ of these two probabilities such that (3.25) 

(k) (k) 

holds. Q\, being a joint distribution, is defined on the product space X)^' x Xj^' whose 
elements are denoted by {q\,q'^)- 

6.1 Definitions and main results 

Recalling that K\(-;x) := i^A(^i,A<=) ^2,A':; 3;) is defined in Definition 3.2 we denote by 

Ao = Ao(gi,Ac, g2,Ac) := {x € l-^Z'^ n A : KA{qi^Ac,q2,Ac;x) > O} (6.1) 
In order to prove Theorem 3.3 we have to find a coupling Q\ so that there is eg such that 

^ QA(eA(x)^) < eg (6.2) 
xeAo 

We define (recall that Bx{R) is the ball of center x and radius R), 

Ai= U 5^(10-2%,^) n A (6.3) 

and we observe that Ai D Aq, dist(Ao, A^) > 10~^°£+,7. 
We denote by 

n = nA = |n(x, s) G N, x G £_,-yZ'^ n A, s G {1, .., S}^ (6.4) 
and in the sequel we will consider only those n such that for all x G i-^^Z'^DA and s G {1, .., 5}, 



n{x, s) ^(fe) 



pW(,) 



Given n and any subset A C A we will call nA the restriction to A of n. 



Given a subset A c A, we call the following metric on 

dA{qA,qA)= d^i^A^q'i) (6.5) 



<i.toY;)= " (6.6) 

1 otherwise 
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We call i?A(A*)/^') the corresponding Wasserstein distance between two measures n and /x' in 



Ra{iJ',IJ'') =inf J dA{qA,qA)dQ{qA,qA) 



= inf ^(^Ancr'^V^Ancr'^O (6-7) 

where the inf runs over all possible joint distributions (couplings) of fi and fi'. 



In Subsection 6.3 we prove the following Theorem. 

Theorem 6.1. Given A union of Na cubes of there is eo = eo(^A) such that for all 

- (k) 

qi,A'^ G '^Ac ; ^ = 1)2, the following holds. 

Given any n' , n" such that n'^^ = =: nAi fAi defined in (6.3)^, the following holds. 

A\A 



Calling Ai = Ai \ S^^ [Ai], for any two configurations q- a\Ai> i = ^-i"^ on '^^^^^ , we denote 
by q.^Ac = q^^^^^^ U ^i,Ac, i = 1, 2. 
Let dG\{q^^\q- ^c,nAi), i = 1,2 be the probabilities dG\{-\qi^A<') , i = 1,2 conditioned to have 
the configuration in equal to q^ i^c and occupation numbers in Ai given by tt-Ai- 
Then for Aq defined in (6.1) 

i?Ao(f^G^(-ki,Af'^Ai),dG^(-k2,A=''^Ai)) <eo (6.8) 

The next result, proved at the end of Subsection 6.6, deals with the Wasserstein distance Ra^ 

of the distributions of the occupation numbers n that in Theorem 6.1 have been set equal to 
each other inside Ai. For these variables the metric dx defined in (6.6) is replaced by 



dx{n,n') 



a n'{x,s) = n"{x,s),'^s 

1 otherwise 



Theorem 6.2. Given A union of Na cubes of there is ei = €i{Na) such that the 

following holds. Let G\{nA\qi,A'') , i = 1,2 be the marginals of dG^{qA\qi,A^) , i = 1,2 on the 
variables ua defined in (6.4). 
Then 

RA,{dGl{nA\qi,A^),dGl{nA\q2,A^)) < ei (6.9) 



In Subsection 6.7 we show that Theorem 3.3 is a consequence of Theorems 6.1 and 6.2. 
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6.2 Two properties of the Wasserstein distance in an abstract setting 

Let be a complete, separable metric space with distance d{uj,uj') and let i?(/xi,/Xo) be the 
corresponding Wasserstein distance between two measures /Lti and hq. Thus 

/io) = inf / d{u),u)')Q{du,dw') (6.10) 
Q J 

where the inf runs over all possible joint distributions of jii and /Uq. 

Theorem 6.3. Let u he a given positive measure on Let h and v he such that for all 
tG[0,l], 

Zt = j e-^^^'^^+^<'^^^iy{du) < oo, (6.11) 

Set 

mt{oj) = Z^^e-^^^'^'^+*'"^'^'^\ Utidu) = mt{uj)iy{du) (6.12) 

Then 

R{lJ-i, fJ-o) < sup Ut(|a;| 1^1) + )Ut(|a;|)/xt(|t;|)) (6.13) 

o<t<i ^ ^ 

where, after fixing arbitrarily an element too G ft, we have called \lo\ = d{io,iOo). 
In particular, 

R{fii,Ho) <2{sup\uj\) (sup|i;(a;)|) (6.14) 



Proof. Let 

m 



(lo) = min{mi(a;), mo(a;)}, C = 1 — J m{u)v{du) 



P{du>du') = {m{u)5^J-^^ + y;{mi{uj) — rn{uj)\[mQ{uj') — m{u>'y\\v{dui)v{du') 
P is a coupling of /xi and /xq and therefore 

P(M1)A*o) < / d{ijO,io')P{duidu)') < I \u)\{\jni{oo) — rn{u])\ + [rnQ{ijj) — rn{(jjy\)i'{doo) 

= / \u!\\mi{u!) — mo{uj)\u{du) 
Jn 

having bounded d{io,Lo') < \lo\ + \uj'\ and integrated over the missing variable. 

d 

(6.13) is then obtained by writing mi{uj) — mo(a;) = / —mt{u). 

Jo 

□ 

The following estimate is taken from [14]: 
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Theorem 6.4. Let A C be a measurable set, fj, a probability on Q and fj,A the probability 
fx conditioned to A. Then 

R{l^, IJ^a) < 2 sup |tc;| ijl{A'') (6.15) 



Proof. Let 



Q{dLJ,du}') = l^^AlJ'{duj)5^{duj') + li^^A''ii{d^)iiA{du') 

where 5^{duj') is the probabihty supported by co. Let / be any bounded, measurable function 
on ri, then 

J f{oj)Q(dLo, du') = f{u)fi{duj) + J^^ f{u)fi{dv) J fiAidoj') = fi{f) 

J f{J)Q{dw,dw') = J f{uj)fiidu;) + fiiA^) J f{u')^iA{du;') 

Hence Q is a coupling and 

R{fx,fiA)< J d(u,u')Q{duj,duj') < J 1i_j^a''{\lo\ + \uj'\)fi{duj)fXAidu!') 
which proves (6.15). □ 

Eventually, we mention the following elementary property: 

Proposition 6.5. Assume that the distance d satisfies m{d) := inf d{io,io') > 0. Then 
for all probability measures /x, u and for all A G ft 

m{d) ■ \fx{A) - iy{A)\ < R{ii,u) (6.16) 

Proof. Without loss of generality we assume n{A) > ^{A). Remarking that la;^u;' > l^eA — 
loj'eA) we get for any coupling Q of /i, i/ 

m{d){n{A) - v{A)) < J d{ui,ui')l^jLu,'G{dco,duj') (6.17) 

and the proposition is proved by taking the infimum over all possible couplings Q. □ 

Remark 6.6. The proposition above states that Wasserstein distances associated to very par- 
ticular distances d are finer than the total variation distance dTvil^-, ^) '■= sup \lJi{A) — i'{A)\. 

In the following, we will use this property for Ra, remarking that m{dA) = 1- 
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6.3 Couplings of multi-canonical measures 



Here we prove Theorem 6.1. Recalling that Ai = Ai \ [^i]) fix two boundary 
conditions = 9i,A\Ai U 9i,A<=) i = 1,2. We have to compare the marginal distributions 

of dG\{q^_^\q- ^c,nAi), i = 1,2 over the configurations in Aq (i.e. well inside Ai). Since the 

probabilities dG'j^{q^^ \q^ ^c, ?^Al)^ i = 1,2 depend only on the restrictions of q. to 52^^- [Ai] 
where n'{x, s) = n"{x, s) the corresponding occupation numbers in the two measures are all 
equal to each other. We will thus study couplings of multi-canonical measures, hence the 
title of the Subsection. 



It is now convenient to label the particles. To this purpose we use a multi-index p = {Cx,s,j), 
where Cx is the cube of T>^^-^ where the particle is; s is its spin and j G {1, ..,n{x, s)} 
distinguishes among the particles in the same cube with same spin. We call JC^^ the set of 
labels 

^Ai ={P= {Cx,s,j),x G Ai,s = l,...,^,^ G {l,..,n{x,s)]] 

Observe that C^_^ is determined by n^^ and we thus have the same labels for the two measures. 
Given p = {Cx, s,j) G we denote by Vp a vector configuration rp = {rj,s) with rj G Cx- 

We then denote by rc^ = {rp,p G -C^^} a vector configuration in Ai. Analogously we define 
fC^c- We then call Hc^ (r£^ k^A^) energy H^^ ^ defined in (3.14) and with uai fixed 
as above. 
Calling 

dfp(r) = IreCa^dr (6.18) 

we define 

—BHr ire \rr ) 

Pc^^{drc^yc^^) = Z{rc^X'e J] ^^{dr) (6.19) 

Remark 6.7. If A is a measurable subset of K\, then Ca denotes all labels {C,s,j) 

with C C A and the marginal o/ Pc^((ir£^|r£^c) over the unlabeled configurations is the 

original multi- canonical measure in A. 

We will thus prove Theorem 6.1 if we can compare 

P' = Pc^^iVc^^) and P" = Pc^^{-\r'i^^) (6.20) 
by evaluating the Wasserstein distance Rao{P',P")- 

We will use the Dobrushin high-temperature techniques which allow to reduce to a comparison 
of the conditional probabilities of a single variable rp. 

Proposition 6.8 (Dobrushin high-temperature theorem). There is c such that the following 
holds. For all po = iCxo,so,jQ), Cxo C Ai, all pi = {Cx-^,si,ji) and all r'p^ and r'p^ 

supRa,{Pc,,{- \r,r'p^),Pc,,{- \r,r;j) < c7'+"-ldist(C.o,c.,)<7-i (6-21) 
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where r = {rp)p^po^p-, 



Proof. The probabilities to compare have the form 
while 

Pr (dr\r r" ) = - e^-^^^'^+^^Wl dr 

yar \L, 'p-, ) ^^^^ ^„ e t J-reC^o «r 

where W-y{r) = — /3V^(r, r^-J and W^(r) = —f3{y^{r,r'-_^) — V^(r, r^-^)} hence 

\W'^{r)\<(3 sup |VV;(r,r')K-<cV+--ldi,t(c.„,c.,)<7-i 

Proposition 6.8 then follows from Theorem 6.3. □ 

Remark 6.9. From the proof above, we see that the r.h.s of (6.21) is actually proportionnal 
to (3^'^'^°^- . In other terms, the effective temperature of the system is of order and 
thus very high indeed. 



Corollary 6.10. With P',P" defined by (6.20), there is eo such that for all 7 small enough 
the following holds: 

RA,{P',P")<eo (6.22) 

Proof. For po andpi as in Proposition 6.8 we call 5(pO)Pi) = c7'^"^""ldist(C2,Q,C2,i)<7-i (which 
is the r.h.s. of (6.21)). Then there is ^ > such that for all 7 small enough the following 
holds: 

i?Ao(j",n < EE E E <^(fO'fi)---'^(^'"'^') 

< g-?dist(Ao,Af) 

The first inequality follows from the Dobrushin high-temperature theorem (Proposition 6.8) 
while the second one is obvious once ^ S{p,p') < 07"" < 1 (which is satisfied for all 7 small 

p'¥=p 

enough) . □ 
In view of Remark 6.7, the Theorem 6.1 is a straightforward consequence of 6.10. □ 
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6.4 Taylor expansion 

In this subsection we consider the marginal of dG^{q\\q/<^c) on the variables pA = ^Z'^ n/\_, 
nA = {'n{x,s),x G i-^^U^ n A, s G {!,.., 5}. By an abuse of notation we denote also the 
marginal with G^(pA|gA<=). 
Recalling (4.2) we get 

Gi(PAl^Ac) = ^iffJ^e-/^^-..^f(PAl.-A=) (6.23) 



Recalling (5.1) we also define 



The following holds: 



Proposition 6.11. For all gi,Ac,92,Ac € X)^c , 

i?Ai(G°(-|gi,Ac),G0(-|g2,Ac)) < RA,{G*{-\qi,Ac),G*{-\q2,Ac)) + 2c7^ (6.25) 
wii/t r 5wen in (4.13). 

Proof. By (4.13) there is c = c{Na) such that 

I^a^'IpaI^ac) - fiPA-, 9Ac)| < cY (6.26) 
By (6.26) and Theorem 6.3, there is a (different) constant c > such that 

Ra, (G°(-|gAc), G*(-|gAc)) < cY (6.27) 
Hence the triangular inequality implies (6.25). □ 

We will bound RAi{G*{-\qi^A''),G*'{-\q2^A'') by using the triangular inequality to replace the 
two measures by their Taylor approximants. 

We first prove the following result true for any P^^+^-measurable region A. 

Theorem 6.12. For any G '^ac , calling fi = G^(-|^a<=); the following holds. 

There are c > and 5 < 1/2 that verifies (6.29) below, so that, calling pA the minimizer of 

/(pa; gA':) 

fi(^{3xeA,3s:\pA{x,s)-pA{x,s)\>ez'^^^'}) < e"^- (6.28) 
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Proof. Denoting simply A := {3x G A, 3s : \pa{x, s) — pa{x, s)\ > we have 

By Theorem 5.6 we have that 

/(pa, ^aO > /(pa, ^aO + ^ (pa - PA, PA - pa) 

oo 



Thus caUing C = ( ^ e'^t" )^'^^ we get 

n=0 



1/(pa;.-ac) ^ exp{-/3^1^^[pA(y,.)-PA(y,.)]^-/3f^?.n 



< e 

n=0 

< e-/3^1/(/5A;9Ac)g-/3f c'(<?+/^-)'' 
We bound the partition function as follows, with < e a small constant to be chosen later: 

"{\PA{x,s)-pA{=:,s)\<€ez''^^^^^xys} 



Ak) 



so that 



M({|PA(a;,s) -PA(a;,s)| > ij^^^^}^ < exp{ 
Remark now that 

r'' (^^yiogiCe-H"/'-') = ay (logT)^ 

with a = {d/2 - 5){l - a_) > 0, 6 = (1 - a_)25 - («+ + aJ)d and c = Ce~^ > 0. Choosing 
5 such that 6 > 0, i.e. 

which is always possible (see (3.3)), we get 7''(log7)'^ — as 7 — 0. The Theorem is now 
proved with < c < /3 '^~'^^ , which is always possible for e small enough. □ 
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We call j the minimizer of /(•; gj,Ac), i = 1,2. We then let 

A<,i = {pAG^^'^:|pA(x,s)-pA,»(x,s)|<r''/'+^Va;,V4, i = 1,2 (6.30) 



- (k) 

Proposition 6.13. For all qi^\c £ Xj^c , i = 1,2, 

RAi{Gl{pA\qi,Ac),Gl{pA\q2,Ac)) < i?Ai (GX(pa|^1,Ac, ^<,l), GA(pA|g2,Ac, ^<,2)) +206^'^-' 

(6.31) 

where G\{pA\qi,A<=, A<^i) i = 1,2 are the probabilities G^{-\qi^A<=) conditioned to A<^i, i = 1,2. 
Proof. (6.31) follows from Theorem 6.12 and Theorem 6.4. □ 



Analogously to (6.3) we define the following subset of A. 

A2= U 5^(10-3%,^) n A (6.32) 

and we observe that A2 D Ai, dist(Ai, A2) > 10~^^i+. We also have 
Lemma 6.14. Let K be as in Theorem 5.1. Then K{x) > for all x G A2. 

Proof. Let x G A2, by definition of k{x), if = B^{1{)-'^^1+^^) n A'^ = then k{x) = 
fh + 1 > 0. Assume then that A^ 7^ 0. By (6.32) and (6.3) there is xq G Aq such that 

\x - xo\ < (1 + 10-^0)10-20£+^^, thus A^ C A^^^ = S^(10-i°£+,^) n A"^ and therefore A^^ ^ 0. 
By definition of Aq wc then have that g^c H Ax = q'Ac D Ax and also that K(xq) = m + 1 > 
with m > 2 where m is given by max |p*'^~'^^(^a=; '"j •s) ~ pi'^^l ^ [Cm+iXm)- Then 

max 1/3^^" ''^-'(^A.c; r, s) - p^''^ < (rn, that implies that k{x) > 0. □ 

r&A^,se{l,..,S} 

Recalling that pi^A is the minimizer of /(s^ijA'^)) i = 1,2, we observe that in general the 
gradient of -Da/ (see (5.6) for notation), evaluated at pi^A does not vanishes in all A. However, 
by Theorem 5.1 and Lemmas 5.4, 6.14 it follows that -DA2/(Pi,A; 9i,Ac) = 0. 
N being defined by Theorem 4.1, we set A2 = A2 U 62^^ "^[^2] and define 



'pi,a{x,s) if X e e-,^Z'^nA2 
Pi,K{x,s) if x G £_,^Z'^ n (A \ A2) 



Pi{x,s) = { ^ ^ ^ ^ ^ r^dr^,.. . (6-33) 



Thus p\ = pi^A in A2 while /9*(x, s) = Pi,a{x, s) for all x G £-Z'^ n A and Vs. We denote by 
p* the common value, thus 

p*{x,s) = pI{x,s) = p*2{x,s), VxG^_^Z'^nA2,Vs (6.34) 
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We also define the matrix Bi \ with entries: 

o / /A fDlf{pi^A;qi^Ac){x,s,x',s') if x,x' e H A2 ,^ 

BiA[x,s,x,s) = <. _ (6.35) 
[Dlf{piX,qi,A^){x,s,x',s') otherwise 

Observe that Bi^\ = L'^/(/?i^a; 9i,A<=)- We denote by B the two matrices restricted to A2 U 
^out ^[^2] which are then equal; their common entries are then 

B{x,s,x',s') = Bi^a{x,s,x',s') = B2,a{x,s,x',s') "ix.x' e ^-^-yZ'' fl (A2), Vs (6.36) 
We define for i = 1 , 2 

ViiPA; QiA") = (^A/(Pi,A; 9i,A=), [pA - P*i]) + ^ {[PA - p*],Bi^A[pA - P*i]) (6-37) 
and the probabilities 

PiiPA) := -^e-f"'- '^^(^^'«-^=)xA< ,(pa), Zi,A = ^'^"^'*'^^^XA< ,(pa) (6.38) 

where XA is the characteristic function of the set A: 
The following holds: 

Proposition 6.15. For all qi^A<^ G -^^vc , ^ = 1, 2, and for all €2 > if ^ is small enough the 
following holds: 

i?Ai(G^(pA|^i,Ac,^<,i),GA(pA|g2,Ac, A<,2) < Ra^Pi, pi-i) + 2c7'^/^ + €2 (6.39) 
Proof. We Taylor expand /(pa; ^i,A=) and we call TZi the third order. 

T^i ■= /(pa; %,Ac) - /(Pi,A; %,A=) - (DAfiPi,A; qi,Ac), [pa - Pi,A]^ 

-^([PA - Pi,A]'Dlf{pi^A; qi,A-)[pA - Pi,A]) (6-40) 
Observe that in A<^i and for a suitable constant ci 

Piim < ci/3£l^|pA(x,5) -p,,a(x,s)|3 < ci£l(^±)V'-'''/' 

and conclude that the right hand side of the above inequality is estimated by cj'^^^ as soon 
as 6 satisfies 



- — 3a_ — 2a+ 



(6.41) 



which is compatible with (6.29), sec (3.3). 

Since Bi^a = -D^/(Pi,a; ^i,a=) and p^ = pi^A, by applying Theorem 6.3 with v = Pi'LlZi and 
h = p£i if {pa; qi,Ac) - 7^l j we get that 

i?Ai(GX(pA|gi,Ac,^<,i),Mi) <C7'^/^ (6.42) 
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Prom Lemma 6.14 and (i) of Theorem 5.1 we get that given any e2 for 7 small enough. 



([pA - P2,a], -D^/(P2,A; g2,Ac)[pA - P2,a]) ^ ([pA - P2]> ^2,a[PA - P2]) 



< 



2 y[PA - P2,a], (L>a/(P2,A; 92,a0 - S2,a) [PA - P2,a]^ 



< 



+ 



([Pl,A - P2,a],S2,a[P1,A - P2,a])^^ 

([pa - P2,a], (£'a/(p2,a; g2,A0 - -Da/(pi,a; g2,A0) [pa - P2,a])^ 



[Pl,A - P2,a], -B2,a[P1,A - P2,a] ) . 



ce 



-10-30(7^+)':' < £2 



By applying Theorem 6.3 with ^; = /3^1[7^2-i([pA-p^], B2,a[pa-P^])] and /i = /3£1(/(pa; g2,Ac)- 
f ) we get that 

i^Ai (GA(PA|g2,Ac, A<,2), M2) < Ct'^/^ + 62 (6.43) 

By using the triangular inequality we then get (6.39). □ 
6.5 Quadratic approximation in continuous variables 



In this subsection we consider the conditional probabilities /Ui(-|Pi,A\A2)' Pi,A\A2 ^ ^<,2> 
i = 1, 2. Since DA2f{Pi,A', ^i,A=) = 0, and recalling (6.34) and (6.36), we have that 



Pi(PA2|Pi,A\A2) •= 



§ ( [PA2 -P*] >'Ba2 [PA2 -P*] ) + ( [PA2 "P*] ,-B[Pi,A\A2 "PI ) 



X%,i(PA2) 



■^j,A2(Pj,A\A2 



(6.44) 



where B/^^ is the matrix B restricted to A2 and where, as usual, ■^i,A2(Pi,A\A2) the sum 
over p^^ of the numerator on the right hand side of (6.44). 

We compare the probabilities Mi("|Pi,A\A2) with measures pi with the same energy but with 
continuous state space. To define these measures we start by setting some notations. 

By convenience we consider the variables nA2 = ^lpA2) thus = {n{x,s),x G ^-Z*^ n 
A2, s G {1, .., S}). Since /Xj, i = 1,2 defined in (6.44) have support on A<^i, the variables 

are such that 

[n(x,s)-a*(x,s)] G {-M,-M+1,...,m}, a* {x, s) = itp* {x, s) (6.45) 

where M is the integer part of i'^'^'^^ (^S as in Theorem 6.12). 
We call ^ = {^{x, s), x e l^U^ n A2, s G {1, .., S}) with 

i{x, s) = r_^l'^ [n(x, s) - a* (x, s)] (6.46) 
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and we denote by Xm = |C • ^{^^ ^) ^ {—M, —M + 1, . . . , M}|. In this new variables the 
boundary conditions become 

= ^:''^'5[ri.,A\A, - «*], n,,A\A, = ^ift,A\A, (6.47) 

By an abuse of notation wc call /Uj(^|.^*) the distribution of the variables ^ under the proba- 
bilities Mi('lft,A\A2) defined in (6.44), thus 

^^^{m = ^e-4i(^.^A,0+(e.^r)] (6.48) 
where Z{$^*) is the sum over ^ G Xm of the numerator. 

We next introduce variables r = {r{x,s),x G A2,s G {l,..,^}) which take values in the 
interval of the real line: 

r{x, s) G £Z'^^^[-M, M + 1] (6.49) 

and we call 

YM = [r: r{x, s) G r'^/^[-M, M + 1], Va; G As, s G {1, .., 5})} (6.50) 
We next define the probabilities measures on Ym as 

dmm) = ^4zsY«"''^'^''''^''^^^''^'^^^^M(r)^^z:, i = 1,2 (6.51) 

where dr = JJ (ir(x, s) and Zm(^j*) is the integral of the numerator. 

x,s 

Proposition 6.16. For all Pi,A\A2 ^ ^<,i' f€.calling (6.47) the following holds: 

i?Ai(/ii(-ieD),/^2(-ia))) < Ri^,{pim).P2m)) + (6.52) 

Proof. Given ^ G Xm we call C{C} = {r : < r{x,s) - £,{x,s) < C'^^^Vx G A2,Vs} we 
define iJ'(||C) as 

^-H'm*),= f e-Pi-J[Ur,B^,r)+(r,^:)]^^ (6.53) 

and the following probabilities rrii on 

e--f^'(iie*) 

"ii(0 = -H'm£*) ' ^ = 1,2 (6.54) 

By continuity there is a point G C(^) such that 

H'mt) = P[l,{ri,BA,r^) + (r^,^*)] (6-55) 
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Therefore 



H'iO-P[l.i^,B^,0 + i^,m < \\V{ir,B^,r)/2 + Om-'^' (6-56) 

where VV'(r) is the vector defined as the gradient of ijj with respect to the variables r{x, s) 
and II ■ II is the norm of the vector •. 
Since H^AaH < c*^ then 

H'iO -^[1{^,Ba,0 + I < c*^-^SiUz'^' < c*57V^(^) (6.57) 

For 7 small , thus by Theorem 6.3 and the triangular inequality we get 

-RAi(m(-|pi,A\A2)>M2(-|P2,A\A2)) < i?Ai (m-i, 771,2) + 207^^/^ (6.58) 

We now observe that at any coupling Q of pi and p2 we can associate a coupling Q* of mi 
and m2 by setting 

Q*{^',a = Q{c{axc{a) 

To prove that Q* is indeed a coupling of mi and m2 we compute for any function ^ on Xm 

Vg, i?Ai(mi,m2) < ^ ciAi(C',?OQ*(e',C") 



Thus 



We next observe that 



(6.59) 



C(«')xC(«") 



^iAi(z:',r")rfQ(z:',r") 



Taking the inf over the coupling Q in the above inequality and using (6.59), we get that 
i?Ai(mi,m2) < Rai{pi{-\Q),P2{-\Q)), thus (6.58) implies (6.52). □ 
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6.6 Gaussian approximation 

We now extend the measures Pi{-\C*) on Ym to a measures Pi, i = 1,2, on the full Euclidean 
space, thus Pi, i = 1,2 are the Gaussian measure defined by the r.h.s. of (6.51) without the 

last characteristic function. 

Thus letting r = {r{x, s) G M"* : x G A2, s G {1, . . . S}), 

dPiizm = ^e-4ifc^A,z:)+(z:,C)] (6.60) 

with Z{Q) the integral of the numerator. 
The following holds: 

Proposition 6.17. There is S* > such that the following holds: 

RA,{pim),P2m)) < RA,{pim),P2m))+'^i'* (6.61) 

Proof. By the Chebischev's inequality, and recalling that VarPj(-|^|') = ||iJA2ll~^) there is c 
such that 

Pi{{\r{x,s)\ > ei}) < dL-y{!f-r\ ^ = 1,2 

By (6.29) there is c^* > such that 

Pi{yh) < E E > ^-}) < (6.62) 

Since pi is equal to the probability Pi conditioned to the set Yui by using Theorem 6.4 and 
the triangular inequality, we get (6.61). □ 

We are thus left with the estimate of -Rai(-Pl, -P2) that we do next. 

Proposition 6.18. There is 63 > such that the following holds: 

RA,(Pim),P2m))<e3 (6.63) 

Proof. We first observe that from the definition of the Wasserstein distance 

i?Ai(Pi(-|6i),P2(-|62)) =infQ(rAi T^rkJ (6-64) 

where is the restriction of r to Ai, namely r^i G 3^Ai '•= {r{x,s) G M.'^,x G Ai,s = 
1, ... 5}. Thus the inf on the r.h.s. of (6.64) can be restricted to all couplings of the marginals 
Pi,Ai on the set ^Ai of the probabilities Pi, i = 1,2. 
Recalling (6.47) we define 

h = B^l^* = rJ'^Bll (i?[n,,A\A2 - «1) > ^ = 1> 2 (6.65) 
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We call the restriction of the vector bi to the set Ai. 

We next call C the matrix with entries Cij = {B/^^)^j^ i = {x,s), j = {x',s'), x,x' G A2, 
s, s' G {!,••• 5"}}, denotes the restriction to Ai of C~^. 

Then remark that marginals of Gaussian variables are Gaussian themselves, so we get: 

(6.66) 

We use that the Wasserstein distance is related to the variational distance via the following 
relation 

2i?Ai(Pl,Ai,P2,Ai) = \\PiAi - P2,Ai\\ (6-67) 

where 

||Pi,Ai - -P2,Ai II := J \ipirAi - &i,Ai) - ■0(rAi - b2,Ai)\drAi (6.68) 



We now prove that 

||A,Ai-P2,aJ| <2||C^^^||||6i,Ai-62,AiI|l2( ^^^''^ (^•^9) 

i=(a;,s),a;ef_Z'*nAi 

To prove (6.69) we interpolate defining M(i) = th\^A\ + (1 — i)^2,An t £ [0)1]- Then, 
shorthanding M = M{t), 

l.h.s. of (6.69) <'^j^j\ {blA^ - ?'2,Ai, C^'(rAi - M)) |V'(rAi - M)drA, dt (6.70) 
Using Cauchy-Schwartz the r.h.s. is bounded by 
<2||C^^^||||6i,Ai-62,Ai||l2 f J [Yl irix,s)-M{x,s)fy^^^irA,-M)drA,dt (6.71) 

hence (6.69). 

To estimate ||6i,Ai — ^2,AiIIl2, we apply Theorem A.l with C = C" = 1,1 the identity 
matrix, and with A = Ba2, observing that Ba.2{x,s,x's') = whenever |x — x'\ > 'j~^N. 
Thus from (A. 10) and (A. 5), using that ft,A\A2 ^ ^<,i) ^ = 1)2, (6.47) and (6.45) we get that 
there are c and c', such that for all x G Ai and since dist(Ai, A2) > lO^^'^^^. 

|&i,Ai(a;,s) - 62,Ai(a;,s)| = | ^ B^l{x,s,y,s')B{r'^^'^ni^f^\A2{y,s) - n2,A\A2{y,s))\ 

s',j/eA\A2 

< ll^ll^i Y e-'^l^-^'lT < c'£ie-^^°"'°^+ 

s',2/eA\A2 

Thus this inequality together with (6.67) and (6.69) implies (6.63). □ 



Proof of Theorem 6.2. Recalling the definition (6.38) of the probabilities ^i, and the 
conditional probabilities defined in (6.44), from Propositions 6.16, 6.17, 6.18 we get that for 
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all Pi,A\A2 e ^<,i> * = 1>2, 

^Ai (/^1(-|/^1,A\A2)>^2(-|P2,A\A2)) < 207^^/^ + 2j^' + €3 = €4 

Thus, there is a coupHng ^(n'^^, n'^2|pi_A\A2) P2,A\A2) c>f the conditional probabiHties /ii(-|Pi,A\A2)) 
z = 1, 2 such that 

Q{n'A, + n'^^\pi,A\^2^P2A\^2) ^ 2e4 (6.72) 

We define for all Pi,A\A2 

^/ / // I A fQ(^'A2'"'A2lPl,A\A2>P2,A\A2) if ft,A\A2 ^ ^<,i, ^ = 1, 2 

Vl'^A2''^A2lPl,A\A2'/'2,A\A2j - S ,^0/ / 1- - nj^O/ // i - - ^ ^1. • 

[(iG^(ra^Jpi_A\A2> «i,A<=)aGA(riA2lP2,A\A2> 92,A=) otherwise 

We then define a coupling Q of the measures fii by letting 

Q("-A2'^42) = Q("'A2'^A2Ipi,a\A2>P2,a\A2)'^Ga(pi,a\A2|9i,aOc^Ga(p2,^^^ (6.73) 
From (6.27), (6.72) and Theorem 6.12 it follows that 

<5«i 7^ "Ai) < 264 + 2cY + 2e-'^- = 65 (6.74) 
Observe that (6.74) implies that 

^Ai(m,M2) < £5 (6.75) 

Then, (6.75), Propositions 6.11, 6.13, 6.15 implies (6.9). □ 

6.7 Proof of Theorem 3.3 

We need to construct a coupling Q\ such that (6.2) holds. 

Recall Ai = Ai \ S^^ [Ai] and that for any two configurations q^^ i = 1,2 on ^^^^ 

we denote by = q^ ^ 'liA"-' ^ = 1)2. Prom Theorem 6.1 we have that, for any nAi, 
there is a coupling QaiI^Ai'^A l^i A'='^2,Af j'^Ai) of the two conditional Gibbs measures 
dG^(g^Jgj^^c,nAi), z = 1,2 such that 

QaM^ ^ ^-'"^ ^ 9a n d'-"^|gi,Af'^2,Aj,«Aj < 2eo (6.76) 

a;e^_,T,Z''nAo 

Given n' andn", we define a coupling = Qai (^Ai^lJ^i.Af ^2,Af' ^''^") of '^^a('I^i,a='^')> 
'^<^a(-|'^2 A<='^^")) by setting 



Qai if '^Ai = ^A 



^1 |dG^(-|g^^c,n')<^^A('l'^2 Af'^") otherwise 
Prom Theorem 6.2 there is a coupling Q* of G°(nA|%,Ac), i = 1,2 such that 

Q*(n'Ai /"aJ <2ei (6.77) 
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Then the final coupHng Qa is defined as follows: 

(^GA(g^^^ J(Zi,A<=, nO<^GA(g^^^Jg2,A'=, n'')<3*(^^'> n") (6.78) 
Thus from (6.76) and (6.77) we get 

<5a (gA n c'i-'^^ ^q'in d^--''^) < eg (6.79) 

To complete the proof of (6.2) we need to show that 

s 

Y Q(9And'-'^^ = g^(nC^-\|p(^--)(g^;x,s)-pW|>Cx(.;,)_i) <e (6.80) 
s=l xeAo 

Since in the set on the l.h.s. of (6.80), q'j^ = q'^, by using (6.27) we have 

QAiq'And'-'-^ = qlndJ-'-\ \p^'-'''\q'j,;x,s)-pi''^\ > Ck(.;x)-i) 

< Gi{\p^'--'\q'j,;x,s)-pi^^ > Ck(.;.)-i;9a) 

+ G° (|/--)(gX;x,.) > CKi.;.)-i;qA) 

< GX(|/--)(c/;;.x,.s) > CKi,.)-i\qAc) 

+ Gli\p^'-'^\q'i;x,s) - pW| > CK(.;x)-i\ql) + 2c7" (6.81) 

Prom Theorem 6.12 and (ii) of Theorem 5.1 it follows that for all x G Aq and for gAc = g^c 
or q'lc), 

GW--Hqa;^,s) - pi'^l > CKi.;.)-i\qAo) < e-^- (6.82) 
which together with (6.81) proves Theorem 3.3. □ 



Part III 

Disagreement percolation 

In this part wc fix t G [0, 1], a bounded P^+'T -measurable region A, /c € {1, . . . , S + 1}; p' and 
p" stand for the measures dG'A('7A, El^^c, T'ac) and dGA{qA,Ii\qAc,T^' \c). They are obtained 
by conditioning measures ly' and v" which could be either DLR measures or Gibbs measures 
dG\' {q\' ,T\q(^A')'') with A' D A. We will first construct a coupling of p' and p" and, with the 
help of such a coupling, we will then define a coupling of ly' and ly" proving that it satisfies 
the requirements of Theorem 3.1. The notation which are most used in this part are reported 
below. 

Main notation and definitions. 

We call 

^ = {q,r) e Xi^^ X Ba (6.83) 
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Given a P(^+'t) measurable subset A of A and ^ = {q,T), we call = (^A^E^) its restriction 
to A. Namely if T = (r(l), . . . T{n)), then 

TAii) = (sp[r(i)] n A, %p[r(i)]nA) , La = (^^(1), ■ ■ • rA(n)) (6.84) 

We will say that we vary ^ in A'^ if we change ^ leaving ^a invariant. 
We denote by f^A the product space, 

nA = {xi'^ xBAf, u; = {C,C')enA (6.85) 

Given a subset A C A and lo = (^, ^') G JIa, we call loa = (■^A, ?a) ^ its restriction to A. 
We call J-\ the (T-algcbra of all Borel sets in JIa and for any I?(^+'t) measurable set A in A 
we call J^A the a-algebra of all Borel sets A such that 1a ('^) does not vary when we change 
u in A'^. 



7 Construction of the coupling 

The target of this section is to construct a "good" coupling Q of fi' and /x". The basic idea 
is to implement the disagreement percolation technique used in van der Berg and Maes, [3], 
Butta ct al., [6], Lebowitz et al,[13]. The first step is to introduce a sequence of random sets 
An, which is done in the next subsection. We will then introduce the notion of "stopping 
sets" and "strong Markov couplings" showing that the sets A„ are indeed stopping sets and, 
using the strong Markov coupling property, we will finally get the desired coupling of /x' and 



7.1 The sequence A„ 

We will define here for each oj = {^',^") € ^Ia a decreasing sequence of P^^+'^^-mcasurable 
sets A„, which are therefore set valued random variables. We set Aq = A and for n > 0, 
define A^+i = A„ \ Sn+i, thus the sequence is defined once we specify the "screening sets" 
E„. Screening sets are defined iteratively with the help of the notion of "good" and "bad 

cubes" . 

After defining in an arbitrary fashion an order among the P^^+'t) cubes of 5q^^^[A], for any 
X){^+,7)-measurable set A c A, we start the definition by calling bad all the cubes of ^out^i^o]- 
We then select among these the first one (according to the pre-definite order) which intersects 
a polymer (i.e. either sp(r') H C 7^ 0, or sp(r") n C 7^ 0), if there is no such cube we then 
take the first cube in ^^^^^[Ao]. Call Ci the cube selected with such a rule. We then define 
^1 ~ '^out^[^i] ^ -^0 and call bad all cubes of Si if Ci intersects a polymer. If not, we say 
that a cube (7 G Ei is good if sp(r') n C = sp(r") n C = and if 

Pi OAo(a;), 6 Ao has been defined in (3.24), (7.1) 
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otherwise C G Ei is called bad. In this way each cube of Si is classified as good or bad and 
therefore all cubes of ^^^^^[Ai] are classified as good or bad. We then select C2 in ^^^^^[Ai] 
in the same way wc had selected Ci in S^^^ [Ao], 112 = n Ai and the cubes of are 

then classified as good or bad by the same rule used for those of Si. By iteration we then 
define a sequence which becomes eventually constant, as it stops changing at A„ if 5out^[A„] 
has no bad cube or if A„ is empty. Since A has N* := |A|/£^ ,^ cubes, A„ is certainly constant 
after N*, but maybe even earlier. In Appendix B we will prove: 

Theorem 7.1. If the sequence {A„} stops at n = N and is non empty, then 

QA n C; [A^] = q'in [Aiv] (7.2) 

and 

sp(r') n [Aiv] = sp(r'0 n 62,^ [An] = % (7.3) 



7.2 Stopping sets 

The random variables A„ are "stopping sets" and the sequence A^j is decreasing, A^_|_i ^ A^j, 
in the following sense. 

• .Fac, a a p(^+.T^) measurable subset of A, is the a algebra of all Borel sets A such that 
1^(0;) does not change if we vary uj in A. 

• A random variable TZ with values in the P(^+>t) measurable subsets of A is called a 
stopping set if for all A, 

{uj^Vl: n{u) = A} G J^Ac (7.4) 

• Two stopping sets IZ' and TZ are such that IZ' ^TZli 

7^'(a;) C ■7^(w), for aU G 

{a; : 7^'(w) = A'} n {u : 7^(u;) = A} G J^a-, for all A' C A 



7.3 Strong Mcirkov couplings 

A coupling Q{duj) of ij! and jj," is called strong Markov in 7?., 7?. a stopping set, if the measure 

dQ{uj) := ^ l{n{oj)=A}d-n-A{'^A\ojA'^)dQ{u}A<^) (7.5) 

AcA 

is also a coupling of n' and /j," for all couplings (i7rA(wA|'^A':) of d/j,' {^a\(,a) , and ^^"(^aI^a)- 
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Theorem 7.2. Given any stopping set IZ, let Q he a coupling of // and fj," which is strong 
Markov in TZ, Then any coupling Q defined by (7.5) is strong Markov in TZ' provided the 
stopping set TZ' is such that TZ' ^ TZ, 



Proof. We have to prove that for any family of couphngs {t: /\{duj /\\uj , IS. C A, u;a<= £ ^^a<=}; 
the probabihty Q{duj) defined as 

dQ{uj) := ^ l{'n'{oj)=A}dv:/\{'^A\i^A'=)dQ{oJA'=) (7.6) 

AcA 

is a coupUng of /i' and fx" . We thus take a function /(^) and wc prove that Q{f) = m'(/)) 
where Q{f), l-^'if), is the expectation of / under Q, respectively fi' . 
Using that TZ' is a stopping set we get 



= XI / Hn'{u,)=A}dQ{u;A-)fi'{f\U'^) 

AcA-^^ 

We now rewrite dQ{uj) by using its definition (7.5) and since TZ' ^TZ we get 

Qif) = X X / l{7e(a;)=A}l{7?,'H=A}f^Q(t^A0 / d'KA{oJA\oJA'^)lJ.'{f\U'') (7.7) 
Observe that (recalling A C A) 

/ d'K^{uA\u^c)tx'{f\U) = [ df,'{U\U^)f^'if\U^,U\A) = /^'(/I^Ac) (7.8) 
We insert (7.8) in (7.7) and we get 

Qif) = X / HTZH=A}dQiu^A^)fi'imAc) = fi'if) 

The Theorem is proved. □ 



7.4 Construction of couplings 

We use the sequence {A„} of decreasing stopping sets (in the order =4) and Theorem 7.2 to 
construct a sequence {Q^} of couplings of fi' and fi", the desired coupling will then be , 
where N* = |A|/4,7- The sequence {Q"} is defined iteratively by setting equal to the 
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product coupling: Q'^ = jj! x jj," which, as it can be easily checked, is strong Markov in Aq. 
Then for any n > we set 

d(3''+^(a;A) = ^ l{A„(a;A)=A}C?7rA(a;A|wA\A,'^A=)c^Q"('^A\A) + l{A„(a;A)=0}'^<3"('^A) C^-^) 

where dQ'^ {uj is the marginal of dQ" over {cja'=} and tta, A 7^ 0, is the coupling of 
dl-'-'{i'/\\i'/\c), and ^^"(CaI^Ac) defined next. We distinguish three cases according to the 
values of WA': = (^a<='^a=)- 

• If WAc is such that either sp(r') n + 0, or sp(r"_) n b'-^^;'{l\\f 0, or both, then 
TTA is the product coupling: <i7rA(CA' ?aI^a=) = <^)"'(CaI^Ac)^m'(^aI^Ac)- 

• If WA- is such that sp(r') n (jJ^t^'IA] = sp(r") n <5j;t''[^] = and g' n Ct'[^] = ^' ^ 

'^out [^] then d7rA(Ck) ^aI'^a^) = V(^aI^Ac)<^(^a -^a)^^A' namely cZtta is the coupling 
supported by the diagonal. 

• Finally let (I>a- be such that sp(r') n b^'^^^ [A] = sp(r") n (5^;^'' [A] = but g' n bl2 [A] / 

n [A]. Call T = U (<5f;t^ n a) , ?7 = A \ T. Let dP{(i^, g^, T' , T") = 

'^/f^'l^^/) r'l?A':)^'^"(9{/' ^''ICa':) h^ the product of the marginal distributions of rf/u'(-|^^c) 
and d\j!\-\^'^c) over Af^^^ x Ba- Let be the coupling defined in Theorem 3.3 and 
letting S = {u;a= : T' n (T U ^^^^^ [T]) = T" n (T U jj^t"' [T]) = 0}, we denote by 1= the 
characteristic function of the set S. 

Then we define 

dTTAlwAlt^AO = IsC^^aOc^Qt (gT)«T|9c/>^A'=>9{/>9A0 ^-P(^C/'^f/'£^''£^") 

+[1 - is(a;A=)]V(9A,r'|eAc)V(9A,r"|eA=) 



By Theorem 7.1 the second case above occurs if and only if all cubes of ^^^^^[A] are good, 
while in the third case there are bad cubes in b^^^ [A] so that is non empty. The proof 

that cubes are good with large probability will be based on Theorem 3.3 and the following 
lemma: 

Lemma 7.3. Suppose A„(a;) = A and that the third case above is verified, namely u^c is 
such that sp(r') n [A] = sp(r") n jj^i^ [A] = and q' fl 5l~l [A] ^ q" n 51,1 [A] ■ Suppose 
also that T' n (T U ^ou^iT]) = T" n (T U S^^+i'[T]) = 0. Let C in S„+i, then C is good if 
WA € @t{x) for all x & C, Qt as in (3.24). 

Proof. The proof follows from the definitions of good cubes and Qt{x) because for all x G C, 
Orix) = Ga(x). □ 
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8 Probability estimates. 



Recall from the beginning of Part III that n' and /i" are obtained by conditioning to the 
configurations outside A the measures v' and v" which are either DLR measures or Gibbs 
measures ^Ga'C^A', r|g(yV')c) with A' 3 A. Thus if Q^* is the coupling of ji' and /x" defined 
in Subsection 7.4, we obtain a coupling P oi v' , v" by writing 

dP{u) = dv'{^'^.)du"{^l.)dQ''\uK\oJ^c), uj = (u;a,cDac), cDac = (Cac,^ac) (8-1) 

We will prove here that there is a constant c such that for all 7 small enough, for any 
p(^+,7)-measurable subset A of A: 

, N _ dist(A,A'^) 

p(^{a;:AAr*(a;)DA}j >l-cie ^+.7 (8.2) 
This proves that ((?^,r') and (g^/,r") agree in A, in the sense of (3.20), with probability 

dist(A,A'=) 

> 1 — cie *+'T from which Theorem 3.1 follows. Indeed if v' and u" arc two DLR 
measures, by the arbitrariness of A and A, (8.2) shows that v' = v" , hence that there is 
a unique DLR measure. If instead v' and v" are two Gibbs measures dGf^i{q\i,V\q\ic) and 
dGA"(gA",r|gA"=), A c A', A c A" then (8.2) yields (3.21). 



8.1 Reduction to a percolation event 

Denote by .4 = A{uj) the union of all bad cubes contained in A and of the cubes in 5^^^ [A] 
with a polymer, namely those cubes C such that C C sp(r), L in L' U L". Since by its 
definition any screening set is connected to a bad cube and since any bad cube in A is 
necessarily contained in a screening set, it follows that if .4 7^ then it is connected to A'^. 
Since the event in (8.2) is bounded by 

{u : Kn*{uj) 3 A}^ C {A{uj) n A 7^ 0} (8.3) 

, it is therefore also bounded by the event that the bad cubes percolate from A to A'^. Hence, 
denoting in the sequel by ^ a connected, P^^+'T^ -measurable subset of A U (^^^^^[A], 

p({Ajv*DAr)< E P{{^ = A}) (8-4) 

We write A = AiL)A2^A3, A; the union of cubes of "type i" . Cubes of type 1 are those with 
a polymer, namely C is type 1 if there is T in F' U F" such that C C sp(F). C is type 2 (also 
called unsuccessful) if C, say in is bad and all cubes of (5q^^^[A„] are without polymers 

(in the above sense). Cubes of type 3 arc the remaining ones, they are therefore in the union 
of ah S„_|_i with Sn+i connected to a type 1 bad cube. Then calling Na = , 

l.h.s. of (8.4) < E 3''\,uS3-A^(n^'^^ = "^^^) ^'-'^ 

^6^+.7^'nA A:A^x,An5'+C [A]^0 
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Since ^3 C (j S'+CIC], 
CeAi 

Na^ < S'^iV^i (8.6) 
Therefore N^^ + A^^^ + S'^N^, > and 

3 

n{A = ^^} C {a, = A2;Na,>^}u{a, = Av,Na,>^^^^^} (8.7) 
We are thus reduced to estimate for any {Ai,A2, A^), 

P{{A2 = A2}), if Na, > ^; P(Mi = ^1}), if Na, > ^^^^^3,^ (8.8) 



8.2 Peierls estimates 

We bound here P{{Ai = ^1}) where Ai is some given set in A U (^^^(''[A]. Thus each cube 
C C ^1 is either contained in sp(r), T G T' or in sp(r), T G T" (or both). Thus 

P{{Ai = Ai}) < 2^^i max max{i/'(sp(r) D B); i/"(sp(r) D B)} (8.9) 

BcAi,Nb>Naj^/2 

where sp(r) = sp(r). Let B = Ci U ■ ■ ■ U Cn, Ci disjoint cubes of V'^^+'f^ then, since u' 

rer 

and v" satisfy the Peierls estimates, 

iy'{sp{r)DB)< J2 iy'{r3ri,...,rn)< e-^p°'^'^-.7(^ri+-+iVrj 

ri,...,r„,sp(ri)DCi ri,...,r„,sp(ri)3Ci 

< g-CpoiC'^l,^JVB/2^ ^ g-CpoiC^^l,^JVr/2^^^ < 2^fie~'^p°i^'^-.T^«/^ (8.10) 
r:sp(r)3C 

for all 7 small enough. Thus 

p{{Ai = Ai}) < 22^Aie-CpoiC^^l,,iV^i/4 (8.11) 



8.3 Probability of unsuccessful cubes 

We will bound here P{{A2 = A2}). Given any n > we define 

A2,nH = M'^) n AnicoT, C'n(w) = N^^^^^^A^ (8-12) 
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9nii^) = XnH • e'^"^'^\ Xn(w) := 1^2.„(a;)=A2nA„(a;)= (8-13) 

e > will be specified later. We are going to prove that for all n, 

£{9n+i) < Sign) < ■ . .< £{go) (8.14) 

where £ is the expectation with respect to P. Since C A = Aq and A2,n*{'^) = -42 (w), 
we then get from (8.14), 

P{{A2 = A2}) < e^^2 (8.15) 
Recalling (8.1), we set P^' = P and for n < N*, 

dP"(a;) = diy'{^'^c)du'{^'ic)dQ''{uA\cjAc), u = (a;A,a)Ac), cDac = (^^c^ac) (8.16) 
calling f " the expectation w.r.t. P". We have £{gn+i) = £'^^^{gn+\)-, hence by (7.9), 



^ ^ N,^n.' I P«(da;AO[l{A„=A,A„+,=A'}Xn(^) 
A A'cA 



£{an+i) 

'' ' ^ .17) 

X / 7rA(da;A I Cac) 

where Sn+i = A \ A'. The last integral is equal to 1 if ^2 H {A \ A'} = 0, while, if this is 
not the case, by (3.25) 

f 7rA(da;A|eAc) <c(e5 + e-^p°'^'^-.7/2), A2 n {A \ A'} ^ (8.18) 

^^2(wA)3>l2n{A\A'} 

We then get from (8.17), 

S{gn+i) < f"(5n)max{l, 1} (8.19) 

We choose 

e=^(c(e, + e-po:C^^^-../^))'" (8.20) 
so that the max on the r.h.s. of (8.19) is 1 which thus proves (8.14) and (8.15). 

8.4 Proof of Theorem 3.1 

As we have shown at the beginning of this Section, Theorem 3.1 follows from (8.2) that we 

prove here. 

Given e as in (8.20), for 7 small enough we bound the r.h.s of (8.11) as 

P{{Al = Ai}) < 22^Aig-CpoiC^€l,^AfAi/4 < ^N^, (g_2l) 

Prom (8.4), (8.5), (8.8), (8.21) and (8.15) we then get 

p({Ajv* D A}'^) < X Yl 3^^2e^^ 

xee+,^ZdnA A:ABx,Ansl+{' [A]^0 

< 2|A| J2 (3e)" 

^ dist(A,A°) 

that implies (8.2). □ 
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Part IV 

Appendices 



A Operators on Euclidean spaces 

For the sake of completeness we recall here some elementary properties of operators on finite 
dimensional Hilbert spaces used in the previous sections. We call H the real Hilbert space of 
vectors u = {u{i)} with scalar product 

{u,v) = ^u{i)v{i) (A.l) 

i 

where i above ranges in a finite index set on which a distance |i — j| is defined (in our 

applications i stands for a pair (x, s), with x G ^Z'^PlA, s G {1, .., 5"}, and either £ = £_^^ or £ = 
7~^/^, A being a fixed D^^^-^^-measurable bounded subset of R'^. Operators on Ti are identified 
to matrices B = B{i,j) by setting Bu{i) = > B{i,j)u{j). We write \u\oo = iaax.\u{i)\, 

j 

2 {Bu,Bu) \Bu\oo 
||S||^ = sup^^ ||-B||oo = sup ' , , ' (A.2) 

Recall that 

||5||oo<max^|S(z,j)|, < max{^ ^ =: \B\ (A.3) 

j j j 

The first inequality in (A.3) is obvious. To prove the second one we write 

^(^i?(i,iMj))'< E \B{i,ji)\\B{z,j2)\l{uUif + uU2f) 

< J2 \Bii,h)\\Bii,h)\uijif <\B\''J2<i)' 



In Theorem A.l below we consider matrices of the form B = C'A~^C", thus including 
{QAQ)~^ (after restricting to QH) and PA{QAQ)~^QA, the matrix considered in (5.37). 
With in mind these two applications we will suppose the diagonal elements of A strictly 
positive and large. 



Theorem A.l. Let B = C'A~^C" with A = D + R, D a diagonal matrix, and suppose there 
are c > 0, c' > and b > such that the following holds (recall the definition of the norm 
\C\ given in (A.3) ). 

\C'\ + \C"\ + \R\<c (A.4) 
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The diagonal elements D{i,i) of D are such that D{i,i) > b for every i. Finally C'{i,j) = 
C"{i,j) = R{i,j) = whenever \i — j\ > c'j~^. Then if b is large enough, 



\B\\ < —, ||i?||oo < max^ j)|e^l''-^'l < (A.5) 



Proof. By (A. 3), \\R\\ < c. On the other hand Iji^ll ^ <b ^ and for b so large that b < 1 
the sum on the r.h.s. of (A.6) below converges and 

CXD 

A-^ = D-^ - D-^RD-^ + D-^RD-^RD-^ = X] ( ~ "-D"^ (A.6) 

n=0 

as seen by multiplying the r.h.s. of (A.6) from the left by A: we then get AD~^ {l—RD~^-\-- ■ ■ ) 
which is equal to 1 after writing AD~^ = 1 + RD~^ and after telescopic cancellations. Thus 
(A.6) holds and 

oo ^ 

p-ii<x:r-'iiiir<^j^ (A.7) 

hence, recalling (A. 3), we get the first inequality in (A.5). We write 

^ |S(i, j)|e^l*-^'l < ^ |C'(i,ii)|e^l*-^il ^ \A-\iui2)\e^^''-''^ ^ \C"{i2,j)\e^^''~^^ 



»2 



< cV"'maxy|A-\zi,Z2)|e^l^^-'^l 
Since ^ \R{ii,i2)\e'^^''~'^^ < e^'c, by (A.6) 



22 



J2\A-\ii,i2)\e^^''-''^ < ^6-"-i[e^'c]^ 



n=0 



hence the second inequality in (A.5). 

□ 



In the next two theorems we consider a matrix Ri with small norm, it represents in our 
applications the matrix PA{QAQ)~^QA which by Theorem A.l has indeed a small norm (if 
b is large). 



Theorem A.2. Let B = A + Ri; suppose A symmetric, {u, Au) > k{u,u) for all u; \\Ri\\ <e 
and K > e > 0. Then B is invertible and 

< -7, K' = K-e (A.8) 
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Suppose further that 

sup^ j)|eTl^-j| < a < oo (A.9) 
' 3 

then 

l^-.(,,.)|<,l + ^)exp{-5^} {A.10) 
Proof. By the integration by parts formula, 

Jo 

Since ||e~"^*|| < e""***, 

||e-^*|| < e-*^* + e-'^* ^ ^ < e-^"-'^' (A.12) 

n=l 

/•oo 

Then / e^^^* is well defined and equal to _B ^; (A. 8) also follows. 
Jo 

Calling ei the vector with components ej(j) = li=j, 

pT poo 

B-\i,j) = (e,,e-^S)+ J (ei,e-^S) (A-13) 

By (A.12), 

/•oo o-k't 

I / (ci, e-^'ej) I < = K - e (A.14) 

Jt ^ 

By a Taylor expansion: 



n=0 ii,..,i„_i 



^^■)I^E^«"^""'' E |S(^,n)|e^l'"'^l---|5(4-i,i)|e^l^-^"-i| (A.15) 
hence using (A.9), 

I / (ei,e-^S-)l< (A-16) 

— ?l 

By choosing r = ^ we then get (A. 10) from (A.14) and (A. 16). 

(X ~\~ K 

□ 



Theorem A. 3. Let B = A + Ri as in Theorem A. 2; call D the diagonal part of A, Rq := 
A — D, R = Rq + Ri and suppose that ||-R||oo < oo. Then 

l|B-'IU<i + ^^l + HU^ (A.17) 

K \ K — e / 
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Proof. Recalling that B = D + R, we use the identity 



B-^ = D-^ - D-^RD-^ + D-^RB-^RD-'^ 



Then 



B-^{i,j) = iei,D-^ej) - {D-^euRD-'ej) +Y,{ei,D-'Rek)iek,B-'eh)ieh,RD-'ej) 

k,h 

SO that 

^ \B-^{i,j)\ < + K-'^\\R\\oo + 



oo 



and (A. 17) follows using (A. 8). 

□ 



B Proof of Theorem 7.1 

In the sequel cubes are always cubes in D(^+'t) and a cube C is called "older" than C if there 
is n such that C C and C C A^. We will prove the theorem as a consequence of the 
following property: 

Property P. Let C be a good cube, x G l-^^.U^ n C , {Ci} the cubes older than C which 
intersect Bx{2'^10~^^£^^^) . If either {Cj} is empty or if all Ci are good, then q'j^ n Cx~'~'^ = 

Before proving Property P, we will use it to prove Theorem 7.1. Suppose that for some N, 

Ajv is non empty and that all cubes in (^^^^^[AAr] are good (thus the sequence A„ stops at 
N). Let C be a cube in ^^^^''[AAr], x € n C and at distance < from A^. Then 

-Ba;(2'^10^"'^*^^-l-.^) n A^ intersects only cubes of (5^^^.^ [Ajv], which are by assumption good; then 
by Property P, n Cx~'''^ — i'a^ Cx^'''\ hence (7.2). (7.3) holds because all cubes of 
^iut [^n] are good. 



We start the proof of Property P by introducing a new function M(.x), x G l-^^'L^. We set 
M{x) = 00 outside A and at all x which are in bad cubes. The definition of M{x) on the 
good cubes is given iteratively in A^. We thus suppose to have already defined M{x) on all 
cubes of A^ and have to define it on S„+i = A'^_^_i \ A^. Let thus C C and x E C. We 

set M{x) = if Ba;(10"^°^+,-y) n A^ = 0, otherwise 

M{x) := 1 + max {M{y)\y G i-^^Z"^ n Bx{W~^°£+,^), y such that (7^^+'^^ C A^} (B.l) 
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To compute the value of M{x), x E C, C C Sn+ij we need to look at all sequences yi,y2, ■■■■ 
such that: \yh — Vh-i] < 10~"^'^^+,7, Cyl^'"^^ is older than Cy^t^^\ ho = x and to know whether 
the cubes Cy^^'"^ are good or bad. In principle the sequence may be arbitrarily long but in 
fact it is not: 

Lemma 1. Let C he a good cube, x G C, then the value of M(x) depends only on whether 
the cubes {Ci} are good or bad, where {Q} is the collection of cubes older than C which 
intersect B^{2'^10-^^e+^^). 

Proof. Since any ball of radius {2'^10~^^ + l)^+,7 intersects at most 2'^ cubes of the partition 
P^^+'t), then any sequence yi,y2, ■■■■ as above consists at most of 2*^ elements. □ 

Since m = 2*^ + 2, then 

either M{x) < m - 2 or M{x) = +oo (B.2) 

We will next prove: 

Lemma 2. Let C be a good cube, x e C, then, if fh — M{x) = h> 0, 

q'a n d'-'-'^ =q'in max {p^'-'-^Xq'j,; x, s) - pW| < ^ (B.3) 

se{i,..,5} 

Proof. The proof is by induction on the "age" of the cubes. We thus suppose that the 
above statements holds for all cubes of A^. Let C be a good cube in then the above 

properties hold by the definition of the function K and of good cubes. □ 

Property P is then an immediate consequence of Lemma 2 and (B.2). 



C Mean field 

In this appendix we prove Theorems 2.2 and 2.3. Our approach is based on the recent works 
[12, 10], having in mind that in [12] the total density was set to 1, the temperature being the 
free parameter, while here we fix the (inverse) temperature /3 = 1, the total density x being 
the free parameter. The two approaches are equivalent, see (2.3). 

To achieve our goal, we will need Lemmas C.l, C.2, C.3, C.4 and C.5 below. The first lemma is 
an essential property relating the total density x to the corresponding constrained minimizer 
in a one-to-one way. The second and third lemmas respectively deal with the first and second 
derivatives of the free energy. They show in particular that the sign of the second derivative 
depends on the roots of some peculiar second degree polynomial. The fourth lemma studies 
the locations of these roots, while the fifth and last lemma gives a general condition for a 
piecewise-convex function to have a common tangent at two different points. 
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The section is organized as follows. We first give some notations and reformulate known 
results, before stating our auxiliary lemmas. Then we prove Theorems 2.2 and 2.3, while the 
proofs of lemmas are deferred to the end of the present section. 



Notations 

For any x G (0, +oo), z G [0, 1], we will denote by the density vector p defined as follows: 

-\ 4^x for i = 2,...,S. ^^-'^ 

Notice that ^ p\^'^^ = x and rewrite (2.4) as follows: 

/"^f(x) = inf {F'^^(p(^'^));0 < z<l}. (C.2) 

Now, remarking that zx = p[^''^^ — P2^''*'', we adapt a result from [10, 12]. Namely, recalling 
Theorem A.l in [12] or section 3 in [10], and comparing (C.3) and (C.4) below with (A. 10) 
and (A. 22) in [12], we know that for any S > 2 there exists a threshold 

xs:=2^ln{S-l) (C.3) 

such that 

• for all X < xs, the function z _p™^(p(^'^)) reaches its minimum at z = 0; 

• for all X > Xs, the function z i— > i?'™^(p(^'^)) reaches its minimum at z = z{x), defined 
as the largest solution of the equation R{z) = x where 

Hi.):=y-±^: (C.4) 



• at X = xs, the function z i— > i*'™^ (p(^'^)) reaches its minimum at z = and at z = 
The statement above means that we have 

First of all, we will see that 

Lemma C.l (Monotony of R and z). The functions R : z ^ '^f^d z : x ^ z(x) are both 

increasing respectively on [zs,l) and [xs,+oo), where zs = §5f- They satisfy the relations 

Ro z = W[a;s,+oo) ^'^^ Z O R = ■ 

Moreover 
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Lemma C.2. 



Lemma C.3. 



lim(/"'*)'(a;) = -oo and lim = +00, 

x—^O X— »+oo 



lim ir'nx) - lim (r')'(x) = u - ^ hs - 1). 



, S-1 1 



(C.6) 
(C.7) 

(C.8) 



'2 ford 



dx"^ 



S-l 
S 



xz{x) I ^^^^ 



r: 



x-R 



z{x) 



(C.9) 



where R^ denotes the roots of the second degree polynomial Pz{X) := — bzX — Cz given 
by 



bz := 



5(S-2) 

5^ 



I — {S-l){l-z)[l+(s-l)z]- 

According to Lemma C.3, the convexity properties of /""""^ will follow from the position of the 
roots of Pz(x) with respect to x. We will actually prove the lemma below 

Lemma C.4 (Roots of Pz). The roots of the polynomial Pz are such that R~ < < Rf and 

• for any S > 60, and for all z G [zs, 1), Rf > R{z); 

• for any 3 < S < 59, there exists a unique Zg G {zs, 1) such that = R{zg). Moreover, 
Rf < R{z) on [zs, z'g) and Rf > R{z) on {zg, 1). 

Eventually, the following fact will be helpful to analyze the convex envelope of /™^: 

Lemma C.5. Let f : {a,b] ^ R and g : [b,c) ^ R be convex functions with continuous 
second derivatives. If f{b) = g{b) and if inf f\x) < g'{b) < f'{b) < snpg'{x), then there 



x<b 



exists a common tangent to their respective graphs Vf,Tg. 



x>b 



We are now ready to prove our theorems. 

Proof of Theorem 2.2. By (C.8), /^'^ is strictly convex. Let us now study the convexity 
of f°"^. Fixing x > xg we remark that Lemma C.l implies z{x) > zs and R{z{x)) = x so 
that Lemma C.4 gives: 

• for all 5, x > > RZ 

' z{x) ' 

• if S" > 60 then R+^. > x; 

• if 3 < 5 < 59, -Ri^-) < X if X < and R'^r^s > x if x > x^, where x^ := R{zg). 
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Therefore, (C.9) shows that if S > 60 then /"'^ is strictly convex on [xs, oo), while if 5 < 59 
then 7°'''^ is strictly concave on [xsjX^] and strictly convex on [xg, +oo). 

Let us analyze the convex envelope of /™^. 

• If 5 > 60, (C.5) and Lemma C.2 show that Lemma C.5 applies to / = f^^^, g = /°''^, 
a = 0, b = xs and c = +oo. 

• If 5 < 59 we first have to deal with the concave part of f"^'^. We introduce the function 
g defined by 



Since {f°^'^)"{x^) = 0, g is convex and has continuous second derivatives. Moreover, 

on [xs,a:^], the graph of 5 is a line located above the graph of /°'''^ (concavity of 
jord^. gjj^gg ^]^g latter intersects the (convex) graph of f^^^, the graph of g and the 
graph of f^^^ intersect at some point with abscisse b G {xs,Xg). Besides, the concavity 
of f°"^ implies g'{b) = {f°"^y{x*g) < {r'^^Yixs), while the convexity of f^^ implies 
if'^ib) > if'^Yixs). Thus Lemma C.2 shows that Lemma C.5 applies to / = /^'^ 
and g defined above. 

In any case, Lemma C.5 implies that there exists a line Ti which is simultaneously tangent 
to the disordered branch of (at some point X- < xs) and to the ordered branch of 
(at some other point x+ > xs)- The function /™^(x) being strictly convex outside [a;_,x_)_], 
the graph of its convex envelope necessarily coincides with Ti (resp. with the graph of /™^) 
inside (resp. outside) Denoting by Ai the slope of Ti, the convex envelope of 

/™f(a;) = /™^(x) — Xix is horizontal on and strictly convex outside this segment of 

minimizers. □ 

Proof of Theorem 2.3. If /? is a minimizer of F = then x = YlsP^ ^ minimizer 

of /{"Ij so that X G {x-, X-^^}. If x = X- < xs, then p = p^^^^^; if x = x^ > xs, then there 
exists k £ {1, . . . , S} such that p = p^'^^ := r^'^ • where Ti^k exchanges the first and the 

k^^ coordinates. Reciprocally, the above 5 + 1 vectors p'^'^^ are all minimizers of F. Moreover, 



thus proving (2.5). 

We now show the second part of Theorem 2.3 dealing with the Hessian of F. Straightforward 
computations show: 





s 



s 



d'^F 



1 

~{* 

Ps 



dpsdps' p=p^ 
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Since p^^^ is a minimizer, := D'^F{p^^^) is semi-definite positive. Actually, L^^^ is definite 
positive, or else the third order corrections in the Taylor-Lagrange formula would contradict 
the extremality of p^: 

qSF 1 
dpsOptdpu Pi 

Taking an orthonormal basis of eigenvectors, the estimate (2.6) holds with k* > the smallest 
eigenvalue of L^^) , L('^+^) . □ 



This section ends with the proofs of Lemma C.l, Lemma C.2, Lemma C.3, Lemma C.4 and 
Lemma C.5 which are stated at the beginning of the section and used in the proofs above. 

Proof of Lemma C.l. We express R'{z) = and show that g is always positive. Recalling 
(C.4) we have: 



S-l 



+ 



1 



l + (5-l)z l-z_ 
1 1 



1 ^ l + {S-l)z 

~^ '^^ ^ ' 

\ — z 



l-z l + {S-l)z 



In 



l + {S-l)z 
l-z 



1 



We now show that g is always positive: 



l-z 1 + {S-I)z 
1 S-l 



In 



1 + {S-I)z 
l-z ' 

S-l 



(l-z)2 [i + (5-l)z]2 1 + {S-I)z l-z 

_ Sz [2{S-l)z-{S-2)] 

- (l-z)2[l + (5-l)z]2 • 

We see immediately that g' > for all z > 2{S-^i) ' that g increases on [§5f ; !)• On this 
subinterval, g is thus minimal at (^fEf j where it takes the value 

"^"^^ -.-2lniS-l)--^ + {S-l), 



^\S-lJ ^' S-l 

which increases with S, vanishes at S" = 2, and is strictly positive for all 5 > 3. From this it 
follows that g is strictly positive on [f^f > 1)) which implies that R is strictly increasing with 
z. Since R goes to +oo when z ^ 1, Lemma C.l is proved. □ 



Proof of Lemma C.2. 

Since p^^'^^ is the vector (^, . . . , |), equations (2.1), (C.5), (C.l) give for all x < xs- 

S{S-1) /X\2 X / X 



r')'ix) = ^x + \n^. (C.IO) 
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Recalling (2.1), (C.5) and (C.l), /°^^(x) = F{x,z{x)) holds for all x > xs, where 



Fix,z) = 



1 5-1 

2 S ' 



x{l — z) x{l — z) 



S 



In- 



x{l + {S-l)z)^ x{l + {S-l)z) 
H ^ — In — — — X. 



(C.ll) 



Using (C.ll) and recalling that (§f = 0) we have for all x > xs- 



S-l 



\x,z{x) 



z{x) 



X{1 - Z') + {S - 1) 



1 - Z 



+ 



1 + {S-I)z 
S 



■^^ x{l + {S-l)z) ^ ^ 



- 1 



-^x + ln-+ln{l-z) + —. 



From (C.4), we know that x > ^ log thus 



(r')'(^)> 



Sz 



1 l0£ 



1 — z S b 



, X xz 



(C.12) 
(C.13) 



(C.14) 



Prom Lemma C.l, z{x) zs ss x ^ xs thus (C.7) follows from (C.IO) - (C.13). Similarly, 
z{x) — 1 as X — oo, thus (C.6) follows by taking limits in (C.IO) and (C.14). □ 



Proof of Lemma C.3. 

First notice that (C.8) follows from (C.IO). Using (C.12) we get: 

in (-)=- + ^(l--^-2-'.) + — .'In^— ^ 

= - + (l - 2^ - 2xz'z) + ^-^z'zx 

X S ' b 

z' /„, S-l 



X 



(^R' + ^[{l-z^)xR'-zx'']^ 



X Z 



where we used = R'{z{x)) (= R' by abusing notations) and zR' = —R{z)+ 

(from (C.4)). This achieves the proof of (C.9). □ 



Proof of Lemma C.4. 
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roots of Pz 

We notice that the discriminant of Pz 

S[S + {3S - A)z] 

is always positive, so that the two distinct roots of Pz are given by 



R± = 
' 2{S - 1) 



1 + {S-I)z 



(C.15) 



For all positive z we have ^^^f^z — negative while i?+ is positive, 

sign of Rf — R{z) 

Wc will actually analyze the sign of Hs{z) := — ^(-2^)]) showing it is strictly 

monotone and thus vanishes at most once. Using (C.4) and (C.15) we get 



Hs{z) 
H'siz) 



Sz 



2(S- l)(l + (5- l)z) 
52 



S-2+ V^iPz)] - log 



1 + {S-I)z 
1 - z 



(C.16) 



S + 2(25 - 3)z + {S- 2){2S - 3)z'^ + (-1 + 2(5 + 2)z + (25 - 3)^2) ^/A(P^ 
2(5 - 1)(1 - z)2(l + (5 - 1)^)2 ^5(S+g^4)^ 



A{z) + B{z)./A(Pz) 



2(5 - 1)(1 - z)\l + (5 - l)z)^./A(P^ 



(C.17) 



In the formula (C.17) above, the denominator as well as the polynomial A(z) in the 
numerator are clearly positive for all z > 0. Since the polynomial B{z) is increasing 
for z > and since B{zs) > -1 + 2(5 - 1)^5 = 25 - 5 > 0, we deduce that Hg{z) is 
always positive for z G [zs,l)- 

• Hs vanishes exactly once <^=^ 5 < 59 

We now check for which values of 5 the function Hg actually vanishes somewhere on 
1). As 2; ^ 1, the leading term in Hs diverges hke (1- 

Z)-V2, so that Hs{z) — > +00. 
Thus Hs will vanish exactly once if and only if Hs{zs) < 0. 



G(5) = Hsizs) 
S 



5(5 - 2)(5 - 2 + ^5(8-115 + 452)) 



21og(5-l) 



G'(5) = 
and 

G'{S) = 



2(5-1)3 

254 - 1053 + 2752 - 405 + 24 - (452 - 135 + 12)^5(8-115 + 452) 
2(5 - 1)V'5(8- 115 + 452) 



25^ - 105^ + 2752 - 405 + 24 = (452 - 135 + 12)^5(8- 115 + 452) 
(25^ - 1053 + 2752 - 405 + 24)2 ^ _ ^ 452) (452 - 135 + 12)2 
4(5 - 2)2(5 - lf{S'^ - 1952 + 485 - 36) = 0. 
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The last bracket reaches a local (negative) maximum at S = 3 — - 1 .4 and a local 
(negative) minimum at S" = ^9-\/2i7 ^ 11.2. Therefore it has exactly one root S*, is 
negative before this root and positive after it. Numerical computations give S* ~ 16.2. 

From this, we know that G is decreasing on [3,5"*] and increasing on [5*, 00). Since 
G(3) < and since G{S) diverges like +^/S as S ^ 00, we get that G has exactly one 
root S > S*, is negative before it and positive after it. Numerical computations show 
5?^ 59.1. □ 



Proof of Lemma C.5. We will use the notation 

K := {ae[b,c);a>b and Tg{a) n T/ / 0} , 
where Tg(a) denotes the tangent to Tg at a. 

Since f{b) = g{b), we have b G K, and K is non-empty. Besides, by continuity of g' , there 
exists 60 G {b,c) such that g'{bo) = f'{b); since f,g are strictly convex, elements of K are 
bounded from above by bo and a* := supK < bo is well defined. 

Now, let an an increasing sequence converging to a*. By definition, Tg(an) intersects Tf, 
and we denote by x„ the abscisse of the intersection point which is the closest to 6, so that 
f'{xn) > g'{oin) > g'{b). We now show that Xn is a bounded decreasing sequence: 

• On {x > Xn\, Tj is above Tf{xn) (convexity of /), which in turn is above Tg{an) 
(definition of x„), and therefore above Tg{an+i) (convexity of g). Thus Tf may not 
intersect Tg{an+i) after abscisse Xn, and Xn+i < x„. 

• By continuity of /', there exists 61 G (a, b) such that f'{bi) = g'{b), thus f'{xn) > g'{b) 
implies Xn > bi (convexity of /). 

Thus G [bi,b] C (a, 5], and by continuity of/, Tg{a*) intersectsTj at {x*,f{x*)). 

In particular, a* G K and f'{x*) > g'{a*). 

If we had f'{x*) > g'{a*) we could apply the implicit function function theorem to '^{a, x) = 
g{a)+g'{a){x — a) — f{x) to deduce that K contains a neighborhood of a*, thus contradicting 
the maximality of a*. Therefore f'{x*) = g'{a*) and Tg{a*) = Tf{x*) is actually tangent to 
Tf. □ 

References 

[1] P. Baffioni, T.Kuna, I Merola, E. Presutti: A liquid vapor phase transition in quantum 
statistical mechanics. Submmitted to Memoirs AMS (2004). 

[2] J. van der Berg: A uniqueness condition for Gibbs measures with application to the two 
dimensional antiferromagnet Commun. Math. Phys. 152(1993), 161-166. 

[3] J. van der Berg, C. Maes: Disagreement percolation in the study of Markov fields Ann. 
Pro6.22(1994), 749-763. 



71 



J. van der Berg, J.E. Steif: Percolation and the hard core lattice model Stochastic 
Processes and Appl. 49(1994), 179-197. 

A. Bovier, M. Zahradnik: The low temperature phase of Kac-Ising models J.Stat. Phys. 
87 (1997), 311-332. 

P. Butta, I. Merola, E. Presutti: On the validity of the van der Waals theory in Ising 
systems with long range interactions Markov Provesses and Related Fields 3 (1977) 63-88 

M. Cassandro, E. Presutti: Phase transitions in Ising systems with long but finite range 
interactions Markov Processes and Related Fields 2(1996) 241-262. 

A. De Masi, I. Merola, E. Presutti, Y. Vignaud: Coexistence of ordered and disordered 
phases in Potts models in the continuum in preparation 

R.L. Dobrushin, S.B. Shlosman: Completely analytical interactions: constructive de- 
scription J.Stat. Phys. 46(5-6)(1987) 983-1014. 

Hans-Otto Georgii, Salvador Miracle-Sole, Jean Ruiz, Valentin Zagrebnov: Mean field 
theory of the Potts Gas J. Phys. A 39 (2006) 9045-9053. 

Hans-Otto Georgii, O. Haggstrom: Phase transition in continuum Potts models. Comm. 
Math. Phys. 181 (1996) 507-528. 

T. Gobron, I. Merola: First order phase transitions in Potts models with finite range 
interactions J. Stat. Phys., 126 (2006). 

J.L. Lebowitz, Mazel, E. Presutti: Liquid vapour phase transitions for systems with 
finite range interactions J. Stat. Phys (1999). 

E. Presutti: Prom Statistical Mechanics towards Continuum Mechanics, course given at 
Max-Plank Institute, Leipzig (1999). 

D. Ruelle: Widom-Rowlinson: Existence of a phase transition in a continuous classical 
system. Phys. Rev. Lett. 27 (1971) 1040-1041. 

M. Zahradnik: A short course on the Pirogov-Sinai theory. Rend. Mat. Appl. 18, 411-486 
(1998). 



72 



